\(\int \frac {x^2}{(a+a \cos (x))^{3/2}} \, dx\) [181]

Optimal result

Integrand size = 14, antiderivative size = 257 \[ \int \frac {x^2}{(a+a \cos (x))^{3/2}} \, dx=-\frac {2 x}{a \sqrt {a+a \cos (x)}}-\frac {i x^2 \arctan \left (e^{\frac {i x}{2}}\right ) \cos \left (\frac {x}{2}\right )}{a \sqrt {a+a \cos (x)}}+\frac {4 \text {arctanh}\left (\sin \left (\frac {x}{2}\right )\right ) \cos \left (\frac {x}{2}\right )}{a \sqrt {a+a \cos (x)}}+\frac {2 i x \cos \left (\frac {x}{2}\right ) \operatorname {PolyLog}\left (2,-i e^{\frac {i x}{2}}\right )}{a \sqrt {a+a \cos (x)}}-\frac {2 i x \cos \left (\frac {x}{2}\right ) \operatorname {PolyLog}\left (2,i e^{\frac {i x}{2}}\right )}{a \sqrt {a+a \cos (x)}}-\frac {4 \cos \left (\frac {x}{2}\right ) \operatorname {PolyLog}\left (3,-i e^{\frac {i x}{2}}\right )}{a \sqrt {a+a \cos (x)}}+\frac {4 \cos \left (\frac {x}{2}\right ) \operatorname {PolyLog}\left (3,i e^{\frac {i x}{2}}\right )}{a \sqrt {a+a \cos (x)}}+\frac {x^2 \tan \left (\frac {x}{2}\right )}{2 a \sqrt {a+a \cos (x)}} \]

[Out]

-2*x/a/(a+a*cos(x))^(1/2)-I*x^2*arctan(exp(1/2*I*x))*cos(1/2*x)/a/(a+a*cos(x))^(1/2)+4*arctanh(sin(1/2*x))*cos
(1/2*x)/a/(a+a*cos(x))^(1/2)+2*I*x*cos(1/2*x)*polylog(2,-I*exp(1/2*I*x))/a/(a+a*cos(x))^(1/2)-2*I*x*cos(1/2*x)
*polylog(2,I*exp(1/2*I*x))/a/(a+a*cos(x))^(1/2)-4*cos(1/2*x)*polylog(3,-I*exp(1/2*I*x))/a/(a+a*cos(x))^(1/2)+4
*cos(1/2*x)*polylog(3,I*exp(1/2*I*x))/a/(a+a*cos(x))^(1/2)+1/2*x^2*tan(1/2*x)/a/(a+a*cos(x))^(1/2)

Rubi [A] (verified)

Time = 0.35 (sec) , antiderivative size = 257, normalized size of antiderivative = 1.00, number of steps used = 10, number of rules used = 7, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.500, Rules used = {3400, 4271, 3855, 4266, 2611, 2320, 6724} \[ \int \frac {x^2}{(a+a \cos (x))^{3/2}} \, dx=-\frac {i x^2 \arctan \left (e^{\frac {i x}{2}}\right ) \cos \left (\frac {x}{2}\right )}{a \sqrt {a \cos (x)+a}}+\frac {4 \cos \left (\frac {x}{2}\right ) \text {arctanh}\left (\sin \left (\frac {x}{2}\right )\right )}{a \sqrt {a \cos (x)+a}}+\frac {2 i x \operatorname {PolyLog}\left (2,-i e^{\frac {i x}{2}}\right ) \cos \left (\frac {x}{2}\right )}{a \sqrt {a \cos (x)+a}}-\frac {2 i x \operatorname {PolyLog}\left (2,i e^{\frac {i x}{2}}\right ) \cos \left (\frac {x}{2}\right )}{a \sqrt {a \cos (x)+a}}-\frac {4 \operatorname {PolyLog}\left (3,-i e^{\frac {i x}{2}}\right ) \cos \left (\frac {x}{2}\right )}{a \sqrt {a \cos (x)+a}}+\frac {4 \operatorname {PolyLog}\left (3,i e^{\frac {i x}{2}}\right ) \cos \left (\frac {x}{2}\right )}{a \sqrt {a \cos (x)+a}}+\frac {x^2 \tan \left (\frac {x}{2}\right )}{2 a \sqrt {a \cos (x)+a}}-\frac {2 x}{a \sqrt {a \cos (x)+a}} \]

[In]

Int[x^2/(a + a*Cos[x])^(3/2),x]

[Out]

(-2*x)/(a*Sqrt[a + a*Cos[x]]) - (I*x^2*ArcTan[E^((I/2)*x)]*Cos[x/2])/(a*Sqrt[a + a*Cos[x]]) + (4*ArcTanh[Sin[x
/2]]*Cos[x/2])/(a*Sqrt[a + a*Cos[x]]) + ((2*I)*x*Cos[x/2]*PolyLog[2, (-I)*E^((I/2)*x)])/(a*Sqrt[a + a*Cos[x]])
- ((2*I)*x*Cos[x/2]*PolyLog[2, I*E^((I/2)*x)])/(a*Sqrt[a + a*Cos[x]]) - (4*Cos[x/2]*PolyLog[3, (-I)*E^((I/2)*
x)])/(a*Sqrt[a + a*Cos[x]]) + (4*Cos[x/2]*PolyLog[3, I*E^((I/2)*x)])/(a*Sqrt[a + a*Cos[x]]) + (x^2*Tan[x/2])/(
2*a*Sqrt[a + a*Cos[x]])

Rule 2320

Int[u_, x_Symbol] :> With[{v = FunctionOfExponential[u, x]}, Dist[v/D[v, x], Subst[Int[FunctionOfExponentialFu
nction[u, x]/x, x], x, v], x]] /; FunctionOfExponentialQ[u, x] &&  !MatchQ[u, (w_)*((a_.)*(v_)^(n_))^(m_) /; F
reeQ[{a, m, n}, x] && IntegerQ[m*n]] &&  !MatchQ[u, E^((c_.)*((a_.) + (b_.)*x))*(F_)[v_] /; FreeQ[{a, b, c}, x
] && InverseFunctionQ[F[x]]]

Rule 2611

Int[Log[1 + (e_.)*((F_)^((c_.)*((a_.) + (b_.)*(x_))))^(n_.)]*((f_.) + (g_.)*(x_))^(m_.), x_Symbol] :> Simp[(-(
f + g*x)^m)*(PolyLog[2, (-e)*(F^(c*(a + b*x)))^n]/(b*c*n*Log[F])), x] + Dist[g*(m/(b*c*n*Log[F])), Int[(f + g*
x)^(m - 1)*PolyLog[2, (-e)*(F^(c*(a + b*x)))^n], x], x] /; FreeQ[{F, a, b, c, e, f, g, n}, x] && GtQ[m, 0]

Rule 3400

Int[((c_.) + (d_.)*(x_))^(m_.)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(n_), x_Symbol] :> Dist[(2*a)^IntPart[n]
*((a + b*Sin[e + f*x])^FracPart[n]/Sin[e/2 + a*(Pi/(4*b)) + f*(x/2)]^(2*FracPart[n])), Int[(c + d*x)^m*Sin[e/2
+ a*(Pi/(4*b)) + f*(x/2)]^(2*n), x], x] /; FreeQ[{a, b, c, d, e, f, m}, x] && EqQ[a^2 - b^2, 0] && IntegerQ[n
+ 1/2] && (GtQ[n, 0] || IGtQ[m, 0])

Rule 3855

Int[csc[(c_.) + (d_.)*(x_)], x_Symbol] :> Simp[-ArcTanh[Cos[c + d*x]]/d, x] /; FreeQ[{c, d}, x]

Rule 4266

Int[csc[(e_.) + Pi*(k_.) + (f_.)*(x_)]*((c_.) + (d_.)*(x_))^(m_.), x_Symbol] :> Simp[-2*(c + d*x)^m*(ArcTanh[E
^(I*k*Pi)*E^(I*(e + f*x))]/f), x] + (-Dist[d*(m/f), Int[(c + d*x)^(m - 1)*Log[1 - E^(I*k*Pi)*E^(I*(e + f*x))],
x], x] + Dist[d*(m/f), Int[(c + d*x)^(m - 1)*Log[1 + E^(I*k*Pi)*E^(I*(e + f*x))], x], x]) /; FreeQ[{c, d, e,
f}, x] && IntegerQ[2*k] && IGtQ[m, 0]

Rule 4271

Int[(csc[(e_.) + (f_.)*(x_)]*(b_.))^(n_)*((c_.) + (d_.)*(x_))^(m_), x_Symbol] :> Simp[(-b^2)*(c + d*x)^m*Cot[e
+ f*x]*((b*Csc[e + f*x])^(n - 2)/(f*(n - 1))), x] + (Dist[b^2*d^2*m*((m - 1)/(f^2*(n - 1)*(n - 2))), Int[(c +
d*x)^(m - 2)*(b*Csc[e + f*x])^(n - 2), x], x] + Dist[b^2*((n - 2)/(n - 1)), Int[(c + d*x)^m*(b*Csc[e + f*x])^
(n - 2), x], x] - Simp[b^2*d*m*(c + d*x)^(m - 1)*((b*Csc[e + f*x])^(n - 2)/(f^2*(n - 1)*(n - 2))), x]) /; Free
Q[{b, c, d, e, f}, x] && GtQ[n, 1] && NeQ[n, 2] && GtQ[m, 1]

Rule 6724

Int[PolyLog[n_, (c_.)*((a_.) + (b_.)*(x_))^(p_.)]/((d_.) + (e_.)*(x_)), x_Symbol] :> Simp[PolyLog[n + 1, c*(a
+ b*x)^p]/(e*p), x] /; FreeQ[{a, b, c, d, e, n, p}, x] && EqQ[b*d, a*e]

Rubi steps \begin{align*} \text {integral}& = \frac {\cos \left (\frac {x}{2}\right ) \int x^2 \sec ^3\left (\frac {x}{2}\right ) \, dx}{2 a \sqrt {a+a \cos (x)}} \\ & = -\frac {2 x}{a \sqrt {a+a \cos (x)}}+\frac {x^2 \tan \left (\frac {x}{2}\right )}{2 a \sqrt {a+a \cos (x)}}+\frac {\cos \left (\frac {x}{2}\right ) \int x^2 \sec \left (\frac {x}{2}\right ) \, dx}{4 a \sqrt {a+a \cos (x)}}+\frac {\left (2 \cos \left (\frac {x}{2}\right )\right ) \int \sec \left (\frac {x}{2}\right ) \, dx}{a \sqrt {a+a \cos (x)}} \\ & = -\frac {2 x}{a \sqrt {a+a \cos (x)}}-\frac {i x^2 \arctan \left (e^{\frac {i x}{2}}\right ) \cos \left (\frac {x}{2}\right )}{a \sqrt {a+a \cos (x)}}+\frac {4 \text {arctanh}\left (\sin \left (\frac {x}{2}\right )\right ) \cos \left (\frac {x}{2}\right )}{a \sqrt {a+a \cos (x)}}+\frac {x^2 \tan \left (\frac {x}{2}\right )}{2 a \sqrt {a+a \cos (x)}}-\frac {\cos \left (\frac {x}{2}\right ) \int x \log \left (1-i e^{\frac {i x}{2}}\right ) \, dx}{a \sqrt {a+a \cos (x)}}+\frac {\cos \left (\frac {x}{2}\right ) \int x \log \left (1+i e^{\frac {i x}{2}}\right ) \, dx}{a \sqrt {a+a \cos (x)}} \\ & = -\frac {2 x}{a \sqrt {a+a \cos (x)}}-\frac {i x^2 \arctan \left (e^{\frac {i x}{2}}\right ) \cos \left (\frac {x}{2}\right )}{a \sqrt {a+a \cos (x)}}+\frac {4 \text {arctanh}\left (\sin \left (\frac {x}{2}\right )\right ) \cos \left (\frac {x}{2}\right )}{a \sqrt {a+a \cos (x)}}+\frac {2 i x \cos \left (\frac {x}{2}\right ) \operatorname {PolyLog}\left (2,-i e^{\frac {i x}{2}}\right )}{a \sqrt {a+a \cos (x)}}-\frac {2 i x \cos \left (\frac {x}{2}\right ) \operatorname {PolyLog}\left (2,i e^{\frac {i x}{2}}\right )}{a \sqrt {a+a \cos (x)}}+\frac {x^2 \tan \left (\frac {x}{2}\right )}{2 a \sqrt {a+a \cos (x)}}-\frac {\left (2 i \cos \left (\frac {x}{2}\right )\right ) \int \operatorname {PolyLog}\left (2,-i e^{\frac {i x}{2}}\right ) \, dx}{a \sqrt {a+a \cos (x)}}+\frac {\left (2 i \cos \left (\frac {x}{2}\right )\right ) \int \operatorname {PolyLog}\left (2,i e^{\frac {i x}{2}}\right ) \, dx}{a \sqrt {a+a \cos (x)}} \\ & = -\frac {2 x}{a \sqrt {a+a \cos (x)}}-\frac {i x^2 \arctan \left (e^{\frac {i x}{2}}\right ) \cos \left (\frac {x}{2}\right )}{a \sqrt {a+a \cos (x)}}+\frac {4 \text {arctanh}\left (\sin \left (\frac {x}{2}\right )\right ) \cos \left (\frac {x}{2}\right )}{a \sqrt {a+a \cos (x)}}+\frac {2 i x \cos \left (\frac {x}{2}\right ) \operatorname {PolyLog}\left (2,-i e^{\frac {i x}{2}}\right )}{a \sqrt {a+a \cos (x)}}-\frac {2 i x \cos \left (\frac {x}{2}\right ) \operatorname {PolyLog}\left (2,i e^{\frac {i x}{2}}\right )}{a \sqrt {a+a \cos (x)}}+\frac {x^2 \tan \left (\frac {x}{2}\right )}{2 a \sqrt {a+a \cos (x)}}-\frac {\left (4 \cos \left (\frac {x}{2}\right )\right ) \text {Subst}\left (\int \frac {\operatorname {PolyLog}(2,-i x)}{x} \, dx,x,e^{\frac {i x}{2}}\right )}{a \sqrt {a+a \cos (x)}}+\frac {\left (4 \cos \left (\frac {x}{2}\right )\right ) \text {Subst}\left (\int \frac {\operatorname {PolyLog}(2,i x)}{x} \, dx,x,e^{\frac {i x}{2}}\right )}{a \sqrt {a+a \cos (x)}} \\ & = -\frac {2 x}{a \sqrt {a+a \cos (x)}}-\frac {i x^2 \arctan \left (e^{\frac {i x}{2}}\right ) \cos \left (\frac {x}{2}\right )}{a \sqrt {a+a \cos (x)}}+\frac {4 \text {arctanh}\left (\sin \left (\frac {x}{2}\right )\right ) \cos \left (\frac {x}{2}\right )}{a \sqrt {a+a \cos (x)}}+\frac {2 i x \cos \left (\frac {x}{2}\right ) \operatorname {PolyLog}\left (2,-i e^{\frac {i x}{2}}\right )}{a \sqrt {a+a \cos (x)}}-\frac {2 i x \cos \left (\frac {x}{2}\right ) \operatorname {PolyLog}\left (2,i e^{\frac {i x}{2}}\right )}{a \sqrt {a+a \cos (x)}}-\frac {4 \cos \left (\frac {x}{2}\right ) \operatorname {PolyLog}\left (3,-i e^{\frac {i x}{2}}\right )}{a \sqrt {a+a \cos (x)}}+\frac {4 \cos \left (\frac {x}{2}\right ) \operatorname {PolyLog}\left (3,i e^{\frac {i x}{2}}\right )}{a \sqrt {a+a \cos (x)}}+\frac {x^2 \tan \left (\frac {x}{2}\right )}{2 a \sqrt {a+a \cos (x)}} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.14 (sec) , antiderivative size = 185, normalized size of antiderivative = 0.72 \[ \int \frac {x^2}{(a+a \cos (x))^{3/2}} \, dx=\frac {\cos \left (\frac {x}{2}\right ) \left (-4 x \cos \left (\frac {x}{2}\right )-2 i x^2 \arctan \left (e^{\frac {i x}{2}}\right ) \cos ^2\left (\frac {x}{2}\right )+8 \text {arctanh}\left (\sin \left (\frac {x}{2}\right )\right ) \cos ^2\left (\frac {x}{2}\right )+4 i x \cos ^2\left (\frac {x}{2}\right ) \operatorname {PolyLog}\left (2,-i e^{\frac {i x}{2}}\right )-4 i x \cos ^2\left (\frac {x}{2}\right ) \operatorname {PolyLog}\left (2,i e^{\frac {i x}{2}}\right )-8 \cos ^2\left (\frac {x}{2}\right ) \operatorname {PolyLog}\left (3,-i e^{\frac {i x}{2}}\right )+8 \cos ^2\left (\frac {x}{2}\right ) \operatorname {PolyLog}\left (3,i e^{\frac {i x}{2}}\right )+x^2 \sin \left (\frac {x}{2}\right )\right )}{(a (1+\cos (x)))^{3/2}} \]

[In]

Integrate[x^2/(a + a*Cos[x])^(3/2),x]

[Out]

(Cos[x/2]*(-4*x*Cos[x/2] - (2*I)*x^2*ArcTan[E^((I/2)*x)]*Cos[x/2]^2 + 8*ArcTanh[Sin[x/2]]*Cos[x/2]^2 + (4*I)*x
*Cos[x/2]^2*PolyLog[2, (-I)*E^((I/2)*x)] - (4*I)*x*Cos[x/2]^2*PolyLog[2, I*E^((I/2)*x)] - 8*Cos[x/2]^2*PolyLog
[3, (-I)*E^((I/2)*x)] + 8*Cos[x/2]^2*PolyLog[3, I*E^((I/2)*x)] + x^2*Sin[x/2]))/(a*(1 + Cos[x]))^(3/2)

Maple [F]

\[\int \frac {x^{2}}{\left (a +\cos \left (x \right ) a \right )^{\frac {3}{2}}}d x\]

[In]

int(x^2/(a+cos(x)*a)^(3/2),x)

[Out]

int(x^2/(a+cos(x)*a)^(3/2),x)

Fricas [F]

\[ \int \frac {x^2}{(a+a \cos (x))^{3/2}} \, dx=\int { \frac {x^{2}}{{\left (a \cos \left (x\right ) + a\right )}^{\frac {3}{2}}} \,d x } \]

[In]

integrate(x^2/(a+a*cos(x))^(3/2),x, algorithm="fricas")

[Out]

integral(sqrt(a*cos(x) + a)*x^2/(a^2*cos(x)^2 + 2*a^2*cos(x) + a^2), x)

Sympy [F]

\[ \int \frac {x^2}{(a+a \cos (x))^{3/2}} \, dx=\int \frac {x^{2}}{\left (a \left (\cos {\left (x \right )} + 1\right )\right )^{\frac {3}{2}}}\, dx \]

[In]

integrate(x**2/(a+a*cos(x))**(3/2),x)

[Out]

Integral(x**2/(a*(cos(x) + 1))**(3/2), x)

Maxima [F]

\[ \int \frac {x^2}{(a+a \cos (x))^{3/2}} \, dx=\int { \frac {x^{2}}{{\left (a \cos \left (x\right ) + a\right )}^{\frac {3}{2}}} \,d x } \]

[In]

integrate(x^2/(a+a*cos(x))^(3/2),x, algorithm="maxima")

[Out]

-1/9*(4*(3*sqrt(2)*x^2*sin(3/2*x) - 4*sqrt(2)*x*cos(3/2*x))*cos(3*x)^3 - 4*(3*sqrt(2)*x^2*cos(3/2*x) + 4*sqrt(
2)*x*sin(3/2*x))*sin(3*x)^3 + 96*sqrt(2)*cos(2*x)^2*sin(3/2*x) + 96*sqrt(2)*sin(2*x)^2*sin(3/2*x) - 4*((12*sqr
t(2)*x*cos(3/2*x) - (9*sqrt(2)*x^2 + 8*sqrt(2))*sin(3/2*x))*cos(2*x) + (12*sqrt(2)*x*cos(x) + (9*sqrt(2)*x^2 -
8*sqrt(2))*sin(x) + 12*sqrt(2)*x)*cos(3/2*x) + (12*sqrt(2)*x*sin(3/2*x) + (9*sqrt(2)*x^2 - 8*sqrt(2))*cos(3/2
*x))*sin(2*x) - (9*sqrt(2)*x^2 - 12*sqrt(2)*x*sin(x) + (9*sqrt(2)*x^2 + 8*sqrt(2))*cos(x))*sin(3/2*x))*cos(3*x
)^2 - 12*((12*sqrt(2)*x*cos(3/2*x) - (9*sqrt(2)*x^2 - 8*sqrt(2))*sin(3/2*x))*cos(3*x) + 3*(12*sqrt(2)*x*cos(3/
2*x) - (9*sqrt(2)*x^2 - 8*sqrt(2))*sin(3/2*x))*cos(2*x) + 3*(12*sqrt(2)*x*cos(x) + (9*sqrt(2)*x^2 - 8*sqrt(2))
*sin(x) + 4*sqrt(2)*x)*cos(3/2*x) + 243*(sqrt(2)*a^2*cos(3*x)^2 + 9*sqrt(2)*a^2*cos(2*x)^2 + 9*sqrt(2)*a^2*cos
(x)^2 + sqrt(2)*a^2*sin(3*x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) + 9*sqrt(2)*a^2*sin
(x)^2 + 6*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2 + 2*(3*sqrt(2)*a^2*cos(2*x) + 3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*c
os(3*x) + 6*(3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(2*x) + 6*(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(x))*sin(
3*x))*integrate(1/9*(x^2*cos(4*x)*cos(3/2*x) + 4*x^2*cos(3*x)*cos(3/2*x) + 6*x^2*cos(2*x)*cos(3/2*x) + x^2*sin
(4*x)*sin(3/2*x) + 4*x^2*sin(3*x)*sin(3/2*x) + 6*x^2*sin(2*x)*sin(3/2*x) + 4*x^2*sin(3/2*x)*sin(x) + (4*x^2*co
s(x) + x^2)*cos(3/2*x))/(a^2*cos(4*x)^2 + 16*a^2*cos(3*x)^2 + 36*a^2*cos(2*x)^2 + 16*a^2*cos(x)^2 + a^2*sin(4*
x)^2 + 16*a^2*sin(3*x)^2 + 36*a^2*sin(2*x)^2 + 48*a^2*sin(2*x)*sin(x) + 16*a^2*sin(x)^2 + 8*a^2*cos(x) + a^2 +
2*(4*a^2*cos(3*x) + 6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(4*x) + 8*(6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*c
os(3*x) + 12*(4*a^2*cos(x) + a^2)*cos(2*x) + 4*(2*a^2*sin(3*x) + 3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(4*x) + 16*
(3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(3*x)), x) - 324*(sqrt(2)*a^2*cos(3*x)^2 + 9*sqrt(2)*a^2*cos(2*x)^2 + 9*sqr
t(2)*a^2*cos(x)^2 + sqrt(2)*a^2*sin(3*x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) + 9*sqr
t(2)*a^2*sin(x)^2 + 6*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2 + 2*(3*sqrt(2)*a^2*cos(2*x) + 3*sqrt(2)*a^2*cos(x) + sq
rt(2)*a^2)*cos(3*x) + 6*(3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(2*x) + 6*(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*
sin(x))*sin(3*x))*integrate(1/9*(x*cos(3/2*x)*sin(4*x) + 4*x*cos(3/2*x)*sin(3*x) + 6*x*cos(3/2*x)*sin(2*x) - x
*cos(4*x)*sin(3/2*x) - 4*x*cos(3*x)*sin(3/2*x) - 6*x*cos(2*x)*sin(3/2*x) + 4*x*cos(3/2*x)*sin(x) - (4*x*cos(x)
+ x)*sin(3/2*x))/(a^2*cos(4*x)^2 + 16*a^2*cos(3*x)^2 + 36*a^2*cos(2*x)^2 + 16*a^2*cos(x)^2 + a^2*sin(4*x)^2 +
16*a^2*sin(3*x)^2 + 36*a^2*sin(2*x)^2 + 48*a^2*sin(2*x)*sin(x) + 16*a^2*sin(x)^2 + 8*a^2*cos(x) + a^2 + 2*(4*
a^2*cos(3*x) + 6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(4*x) + 8*(6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(3*x
) + 12*(4*a^2*cos(x) + a^2)*cos(2*x) + 4*(2*a^2*sin(3*x) + 3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(4*x) + 16*(3*a^2
*sin(2*x) + 2*a^2*sin(x))*sin(3*x)), x) + (12*sqrt(2)*x*sin(3/2*x) + (9*sqrt(2)*x^2 - 8*sqrt(2))*cos(3/2*x))*s
in(3*x) + 3*(12*sqrt(2)*x*sin(3/2*x) + (9*sqrt(2)*x^2 - 8*sqrt(2))*cos(3/2*x))*sin(2*x) - (9*sqrt(2)*x^2 - 36*
sqrt(2)*x*sin(x) + 3*(9*sqrt(2)*x^2 - 8*sqrt(2))*cos(x) - 8*sqrt(2))*sin(3/2*x))*cos(4/3*arctan2(sin(3/2*x), c
os(3/2*x)))^2 - 12*((12*sqrt(2)*x*cos(3/2*x) - (9*sqrt(2)*x^2 - 8*sqrt(2))*sin(3/2*x))*cos(3*x) + 3*(12*sqrt(2
)*x*cos(3/2*x) - (9*sqrt(2)*x^2 - 8*sqrt(2))*sin(3/2*x))*cos(2*x) + 3*(12*sqrt(2)*x*cos(x) + (9*sqrt(2)*x^2 -
8*sqrt(2))*sin(x) + 4*sqrt(2)*x)*cos(3/2*x) + 243*(sqrt(2)*a^2*cos(3*x)^2 + 9*sqrt(2)*a^2*cos(2*x)^2 + 9*sqrt(
2)*a^2*cos(x)^2 + sqrt(2)*a^2*sin(3*x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) + 9*sqrt(
2)*a^2*sin(x)^2 + 6*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2 + 2*(3*sqrt(2)*a^2*cos(2*x) + 3*sqrt(2)*a^2*cos(x) + sqrt
(2)*a^2)*cos(3*x) + 6*(3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(2*x) + 6*(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*si
n(x))*sin(3*x))*integrate(1/9*(x^2*cos(4*x)*cos(3/2*x) + 4*x^2*cos(3*x)*cos(3/2*x) + 6*x^2*cos(2*x)*cos(3/2*x)
+ x^2*sin(4*x)*sin(3/2*x) + 4*x^2*sin(3*x)*sin(3/2*x) + 6*x^2*sin(2*x)*sin(3/2*x) + 4*x^2*sin(3/2*x)*sin(x) +
(4*x^2*cos(x) + x^2)*cos(3/2*x))/(a^2*cos(4*x)^2 + 16*a^2*cos(3*x)^2 + 36*a^2*cos(2*x)^2 + 16*a^2*cos(x)^2 +
a^2*sin(4*x)^2 + 16*a^2*sin(3*x)^2 + 36*a^2*sin(2*x)^2 + 48*a^2*sin(2*x)*sin(x) + 16*a^2*sin(x)^2 + 8*a^2*cos(
x) + a^2 + 2*(4*a^2*cos(3*x) + 6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(4*x) + 8*(6*a^2*cos(2*x) + 4*a^2*cos(x
) + a^2)*cos(3*x) + 12*(4*a^2*cos(x) + a^2)*cos(2*x) + 4*(2*a^2*sin(3*x) + 3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(
4*x) + 16*(3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(3*x)), x) - 324*(sqrt(2)*a^2*cos(3*x)^2 + 9*sqrt(2)*a^2*cos(2*x)
^2 + 9*sqrt(2)*a^2*cos(x)^2 + sqrt(2)*a^2*sin(3*x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(
x) + 9*sqrt(2)*a^2*sin(x)^2 + 6*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2 + 2*(3*sqrt(2)*a^2*cos(2*x) + 3*sqrt(2)*a^2*c
os(x) + sqrt(2)*a^2)*cos(3*x) + 6*(3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(2*x) + 6*(sqrt(2)*a^2*sin(2*x) + sq
rt(2)*a^2*sin(x))*sin(3*x))*integrate(1/9*(x*cos(3/2*x)*sin(4*x) + 4*x*cos(3/2*x)*sin(3*x) + 6*x*cos(3/2*x)*si
n(2*x) - x*cos(4*x)*sin(3/2*x) - 4*x*cos(3*x)*sin(3/2*x) - 6*x*cos(2*x)*sin(3/2*x) + 4*x*cos(3/2*x)*sin(x) - (
4*x*cos(x) + x)*sin(3/2*x))/(a^2*cos(4*x)^2 + 16*a^2*cos(3*x)^2 + 36*a^2*cos(2*x)^2 + 16*a^2*cos(x)^2 + a^2*si
n(4*x)^2 + 16*a^2*sin(3*x)^2 + 36*a^2*sin(2*x)^2 + 48*a^2*sin(2*x)*sin(x) + 16*a^2*sin(x)^2 + 8*a^2*cos(x) + a
^2 + 2*(4*a^2*cos(3*x) + 6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(4*x) + 8*(6*a^2*cos(2*x) + 4*a^2*cos(x) + a^
2)*cos(3*x) + 12*(4*a^2*cos(x) + a^2)*cos(2*x) + 4*(2*a^2*sin(3*x) + 3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(4*x) +
16*(3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(3*x)), x) + (12*sqrt(2)*x*sin(3/2*x) + (9*sqrt(2)*x^2 - 8*sqrt(2))*cos
(3/2*x))*sin(3*x) + 3*(12*sqrt(2)*x*sin(3/2*x) + (9*sqrt(2)*x^2 - 8*sqrt(2))*cos(3/2*x))*sin(2*x) - (9*sqrt(2)
*x^2 - 36*sqrt(2)*x*sin(x) + 3*(9*sqrt(2)*x^2 - 8*sqrt(2))*cos(x) - 8*sqrt(2))*sin(3/2*x))*cos(2/3*arctan2(sin
(3/2*x), cos(3/2*x)))^2 + 4*((3*sqrt(2)*x^2*sin(3/2*x) - 4*sqrt(2)*x*cos(3/2*x))*cos(3*x) - (12*sqrt(2)*x*cos(
3/2*x) - (9*sqrt(2)*x^2 - 8*sqrt(2))*sin(3/2*x))*cos(2*x) - (12*sqrt(2)*x*cos(x) + (9*sqrt(2)*x^2 + 8*sqrt(2))
*sin(x) + 4*sqrt(2)*x)*cos(3/2*x) - (12*sqrt(2)*x*sin(3/2*x) + (9*sqrt(2)*x^2 + 8*sqrt(2))*cos(3/2*x))*sin(2*x
) + (3*sqrt(2)*x^2 - 12*sqrt(2)*x*sin(x) + (9*sqrt(2)*x^2 - 8*sqrt(2))*cos(x))*sin(3/2*x))*sin(3*x)^2 - 12*((1
2*sqrt(2)*x*cos(3/2*x) - (9*sqrt(2)*x^2 - 8*sqrt(2))*sin(3/2*x))*cos(3*x) + 3*(12*sqrt(2)*x*cos(3/2*x) - (9*sq
rt(2)*x^2 - 8*sqrt(2))*sin(3/2*x))*cos(2*x) + 3*(12*sqrt(2)*x*cos(x) + (9*sqrt(2)*x^2 - 8*sqrt(2))*sin(x) + 4*
sqrt(2)*x)*cos(3/2*x) + 243*(sqrt(2)*a^2*cos(3*x)^2 + 9*sqrt(2)*a^2*cos(2*x)^2 + 9*sqrt(2)*a^2*cos(x)^2 + sqrt
(2)*a^2*sin(3*x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) + 9*sqrt(2)*a^2*sin(x)^2 + 6*sq
rt(2)*a^2*cos(x) + sqrt(2)*a^2 + 2*(3*sqrt(2)*a^2*cos(2*x) + 3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(3*x) + 6*
(3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(2*x) + 6*(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(x))*sin(3*x))*integr
ate(1/9*(x^2*cos(4*x)*cos(3/2*x) + 4*x^2*cos(3*x)*cos(3/2*x) + 6*x^2*cos(2*x)*cos(3/2*x) + x^2*sin(4*x)*sin(3/
2*x) + 4*x^2*sin(3*x)*sin(3/2*x) + 6*x^2*sin(2*x)*sin(3/2*x) + 4*x^2*sin(3/2*x)*sin(x) + (4*x^2*cos(x) + x^2)*
cos(3/2*x))/(a^2*cos(4*x)^2 + 16*a^2*cos(3*x)^2 + 36*a^2*cos(2*x)^2 + 16*a^2*cos(x)^2 + a^2*sin(4*x)^2 + 16*a^
2*sin(3*x)^2 + 36*a^2*sin(2*x)^2 + 48*a^2*sin(2*x)*sin(x) + 16*a^2*sin(x)^2 + 8*a^2*cos(x) + a^2 + 2*(4*a^2*co
s(3*x) + 6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(4*x) + 8*(6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(3*x) + 12
*(4*a^2*cos(x) + a^2)*cos(2*x) + 4*(2*a^2*sin(3*x) + 3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(4*x) + 16*(3*a^2*sin(2
*x) + 2*a^2*sin(x))*sin(3*x)), x) - 324*(sqrt(2)*a^2*cos(3*x)^2 + 9*sqrt(2)*a^2*cos(2*x)^2 + 9*sqrt(2)*a^2*cos
(x)^2 + sqrt(2)*a^2*sin(3*x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) + 9*sqrt(2)*a^2*sin
(x)^2 + 6*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2 + 2*(3*sqrt(2)*a^2*cos(2*x) + 3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*c
os(3*x) + 6*(3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(2*x) + 6*(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(x))*sin(
3*x))*integrate(1/9*(x*cos(3/2*x)*sin(4*x) + 4*x*cos(3/2*x)*sin(3*x) + 6*x*cos(3/2*x)*sin(2*x) - x*cos(4*x)*si
n(3/2*x) - 4*x*cos(3*x)*sin(3/2*x) - 6*x*cos(2*x)*sin(3/2*x) + 4*x*cos(3/2*x)*sin(x) - (4*x*cos(x) + x)*sin(3/
2*x))/(a^2*cos(4*x)^2 + 16*a^2*cos(3*x)^2 + 36*a^2*cos(2*x)^2 + 16*a^2*cos(x)^2 + a^2*sin(4*x)^2 + 16*a^2*sin(
3*x)^2 + 36*a^2*sin(2*x)^2 + 48*a^2*sin(2*x)*sin(x) + 16*a^2*sin(x)^2 + 8*a^2*cos(x) + a^2 + 2*(4*a^2*cos(3*x)
+ 6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(4*x) + 8*(6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(3*x) + 12*(4*a^
2*cos(x) + a^2)*cos(2*x) + 4*(2*a^2*sin(3*x) + 3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(4*x) + 16*(3*a^2*sin(2*x) +
2*a^2*sin(x))*sin(3*x)), x) + (12*sqrt(2)*x*sin(3/2*x) + (9*sqrt(2)*x^2 - 8*sqrt(2))*cos(3/2*x))*sin(3*x) + 3*
(12*sqrt(2)*x*sin(3/2*x) + (9*sqrt(2)*x^2 - 8*sqrt(2))*cos(3/2*x))*sin(2*x) - (9*sqrt(2)*x^2 - 36*sqrt(2)*x*si
n(x) + 3*(9*sqrt(2)*x^2 - 8*sqrt(2))*cos(x) - 8*sqrt(2))*sin(3/2*x))*sin(4/3*arctan2(sin(3/2*x), cos(3/2*x)))^
2 - 12*((12*sqrt(2)*x*cos(3/2*x) - (9*sqrt(2)*x^2 - 8*sqrt(2))*sin(3/2*x))*cos(3*x) + 3*(12*sqrt(2)*x*cos(3/2*
x) - (9*sqrt(2)*x^2 - 8*sqrt(2))*sin(3/2*x))*cos(2*x) + 3*(12*sqrt(2)*x*cos(x) + (9*sqrt(2)*x^2 - 8*sqrt(2))*s
in(x) + 4*sqrt(2)*x)*cos(3/2*x) + 243*(sqrt(2)*a^2*cos(3*x)^2 + 9*sqrt(2)*a^2*cos(2*x)^2 + 9*sqrt(2)*a^2*cos(x
)^2 + sqrt(2)*a^2*sin(3*x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) + 9*sqrt(2)*a^2*sin(x
)^2 + 6*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2 + 2*(3*sqrt(2)*a^2*cos(2*x) + 3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos
(3*x) + 6*(3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(2*x) + 6*(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(x))*sin(3*
x))*integrate(1/9*(x^2*cos(4*x)*cos(3/2*x) + 4*x^2*cos(3*x)*cos(3/2*x) + 6*x^2*cos(2*x)*cos(3/2*x) + x^2*sin(4
*x)*sin(3/2*x) + 4*x^2*sin(3*x)*sin(3/2*x) + 6*x^2*sin(2*x)*sin(3/2*x) + 4*x^2*sin(3/2*x)*sin(x) + (4*x^2*cos(
x) + x^2)*cos(3/2*x))/(a^2*cos(4*x)^2 + 16*a^2*cos(3*x)^2 + 36*a^2*cos(2*x)^2 + 16*a^2*cos(x)^2 + a^2*sin(4*x)
^2 + 16*a^2*sin(3*x)^2 + 36*a^2*sin(2*x)^2 + 48*a^2*sin(2*x)*sin(x) + 16*a^2*sin(x)^2 + 8*a^2*cos(x) + a^2 + 2
*(4*a^2*cos(3*x) + 6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(4*x) + 8*(6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos
(3*x) + 12*(4*a^2*cos(x) + a^2)*cos(2*x) + 4*(2*a^2*sin(3*x) + 3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(4*x) + 16*(3
*a^2*sin(2*x) + 2*a^2*sin(x))*sin(3*x)), x) - 324*(sqrt(2)*a^2*cos(3*x)^2 + 9*sqrt(2)*a^2*cos(2*x)^2 + 9*sqrt(
2)*a^2*cos(x)^2 + sqrt(2)*a^2*sin(3*x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) + 9*sqrt(
2)*a^2*sin(x)^2 + 6*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2 + 2*(3*sqrt(2)*a^2*cos(2*x) + 3*sqrt(2)*a^2*cos(x) + sqrt
(2)*a^2)*cos(3*x) + 6*(3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(2*x) + 6*(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*si
n(x))*sin(3*x))*integrate(1/9*(x*cos(3/2*x)*sin(4*x) + 4*x*cos(3/2*x)*sin(3*x) + 6*x*cos(3/2*x)*sin(2*x) - x*c
os(4*x)*sin(3/2*x) - 4*x*cos(3*x)*sin(3/2*x) - 6*x*cos(2*x)*sin(3/2*x) + 4*x*cos(3/2*x)*sin(x) - (4*x*cos(x) +
x)*sin(3/2*x))/(a^2*cos(4*x)^2 + 16*a^2*cos(3*x)^2 + 36*a^2*cos(2*x)^2 + 16*a^2*cos(x)^2 + a^2*sin(4*x)^2 + 1
6*a^2*sin(3*x)^2 + 36*a^2*sin(2*x)^2 + 48*a^2*sin(2*x)*sin(x) + 16*a^2*sin(x)^2 + 8*a^2*cos(x) + a^2 + 2*(4*a^
2*cos(3*x) + 6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(4*x) + 8*(6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(3*x)
+ 12*(4*a^2*cos(x) + a^2)*cos(2*x) + 4*(2*a^2*sin(3*x) + 3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(4*x) + 16*(3*a^2*s
in(2*x) + 2*a^2*sin(x))*sin(3*x)), x) + (12*sqrt(2)*x*sin(3/2*x) + (9*sqrt(2)*x^2 - 8*sqrt(2))*cos(3/2*x))*sin
(3*x) + 3*(12*sqrt(2)*x*sin(3/2*x) + (9*sqrt(2)*x^2 - 8*sqrt(2))*cos(3/2*x))*sin(2*x) - (9*sqrt(2)*x^2 - 36*sq
rt(2)*x*sin(x) + 3*(9*sqrt(2)*x^2 - 8*sqrt(2))*cos(x) - 8*sqrt(2))*sin(3/2*x))*sin(2/3*arctan2(sin(3/2*x), cos
(3/2*x)))^2 + 4*(24*sqrt(2)*cos(2*x)^2*sin(3/2*x) + 24*sqrt(2)*sin(2*x)^2*sin(3/2*x) - 2*(12*sqrt(2)*x*cos(3/2
*x) - (9*sqrt(2)*x^2 + 24*sqrt(2)*cos(x) + 8*sqrt(2))*sin(3/2*x))*cos(2*x) - 2*(12*sqrt(2)*x*cos(x) + (9*sqrt(
2)*x^2 - 8*sqrt(2))*sin(x) + 6*sqrt(2)*x)*cos(3/2*x) - 2*((9*sqrt(2)*x^2 - 8*sqrt(2))*cos(3/2*x) + 12*(sqrt(2)
*x - 2*sqrt(2)*sin(x))*sin(3/2*x))*sin(2*x) + (9*sqrt(2)*x^2 + 24*sqrt(2)*cos(x)^2 - 24*sqrt(2)*x*sin(x) + 24*
sqrt(2)*sin(x)^2 + 2*(9*sqrt(2)*x^2 + 8*sqrt(2))*cos(x))*sin(3/2*x))*cos(3*x) - 4*(12*sqrt(2)*x*cos(3/2*x) - (
9*sqrt(2)*x^2 + 48*sqrt(2)*cos(x) + 8*sqrt(2))*sin(3/2*x))*cos(2*x) - 4*(12*sqrt(2)*x*cos(x) + (9*sqrt(2)*x^2
- 8*sqrt(2))*sin(x) + 4*sqrt(2)*x)*cos(3/2*x) - 4*(sqrt(2)*sin(3*x)^3 + 6*(sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*
sin(3*x)^2 + (sqrt(2)*cos(3*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(
2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*cos(x) + sqrt(2))*cos(3*x) + 6*(3*
sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*sqrt(2)*cos(x) + sqrt(2))*sin(3*x) + 3*(sqrt(2)*cos(3*x)^2 + 9*sqrt(2)*
cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + sqrt(2)*sin(3*x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*s
qrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*cos(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2)
)*cos(2*x) + 6*(sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*sin(3*x) + 6*sqrt(2)*cos(x) + sqrt(2))*sin(4/3*arctan2(sin(
3/2*x), cos(3/2*x))) + 3*(sqrt(2)*cos(3*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + sqrt(2)*sin(3*x)^2
+ 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*c
os(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*(sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*si
n(3*x) + 6*sqrt(2)*cos(x) + sqrt(2))*sin(2/3*arctan2(sin(3/2*x), cos(3/2*x))))*cos(5/3*arctan2(sin(3/2*x), cos
(3/2*x))) + 4*(72*sqrt(2)*cos(2*x)^2*sin(3/2*x) + 8*sqrt(2)*sin(3*x)^2*sin(3/2*x) + 72*sqrt(2)*sin(2*x)^2*sin(
3/2*x) - 2*(12*sqrt(2)*x*cos(3/2*x) - (9*sqrt(2)*x^2 - 4*sqrt(2))*sin(3/2*x))*cos(3*x)^2 + 2*(9*(3*sqrt(2)*x^2
*sin(3/2*x) - 4*sqrt(2)*x*cos(3/2*x))*cos(2*x) - 3*(12*sqrt(2)*x*cos(x) + (9*sqrt(2)*x^2 - 8*sqrt(2))*sin(x) +
8*sqrt(2)*x)*cos(3/2*x) - 3*(12*sqrt(2)*x*sin(3/2*x) + (9*sqrt(2)*x^2 - 8*sqrt(2))*cos(3/2*x))*sin(2*x) + (27
*sqrt(2)*x^2*cos(x) + 18*sqrt(2)*x^2 - 36*sqrt(2)*x*sin(x) - 8*sqrt(2))*sin(3/2*x))*cos(3*x) - 18*(4*sqrt(2)*x
*cos(3/2*x) - (3*sqrt(2)*x^2 + 8*sqrt(2)*cos(x))*sin(3/2*x))*cos(2*x) - 6*(12*sqrt(2)*x*cos(x) + (9*sqrt(2)*x^
2 - 8*sqrt(2))*sin(x) + 4*sqrt(2)*x)*cos(3/2*x) - 6*((12*sqrt(2)*x*cos(3/2*x) - (9*sqrt(2)*x^2 - 8*sqrt(2))*si
n(3/2*x))*cos(3*x) + 3*(12*sqrt(2)*x*cos(3/2*x) - (9*sqrt(2)*x^2 - 8*sqrt(2))*sin(3/2*x))*cos(2*x) + 3*(12*sqr
t(2)*x*cos(x) + (9*sqrt(2)*x^2 - 8*sqrt(2))*sin(x) + 4*sqrt(2)*x)*cos(3/2*x) + 243*(sqrt(2)*a^2*cos(3*x)^2 + 9
*sqrt(2)*a^2*cos(2*x)^2 + 9*sqrt(2)*a^2*cos(x)^2 + sqrt(2)*a^2*sin(3*x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt
(2)*a^2*sin(2*x)*sin(x) + 9*sqrt(2)*a^2*sin(x)^2 + 6*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2 + 2*(3*sqrt(2)*a^2*cos(2
*x) + 3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(3*x) + 6*(3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(2*x) + 6*(sqrt
(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(x))*sin(3*x))*integrate(1/9*(x^2*cos(4*x)*cos(3/2*x) + 4*x^2*cos(3*x)*cos(3
/2*x) + 6*x^2*cos(2*x)*cos(3/2*x) + x^2*sin(4*x)*sin(3/2*x) + 4*x^2*sin(3*x)*sin(3/2*x) + 6*x^2*sin(2*x)*sin(3
/2*x) + 4*x^2*sin(3/2*x)*sin(x) + (4*x^2*cos(x) + x^2)*cos(3/2*x))/(a^2*cos(4*x)^2 + 16*a^2*cos(3*x)^2 + 36*a^
2*cos(2*x)^2 + 16*a^2*cos(x)^2 + a^2*sin(4*x)^2 + 16*a^2*sin(3*x)^2 + 36*a^2*sin(2*x)^2 + 48*a^2*sin(2*x)*sin(
x) + 16*a^2*sin(x)^2 + 8*a^2*cos(x) + a^2 + 2*(4*a^2*cos(3*x) + 6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(4*x)
+ 8*(6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(3*x) + 12*(4*a^2*cos(x) + a^2)*cos(2*x) + 4*(2*a^2*sin(3*x) + 3*
a^2*sin(2*x) + 2*a^2*sin(x))*sin(4*x) + 16*(3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(3*x)), x) - 324*(sqrt(2)*a^2*co
s(3*x)^2 + 9*sqrt(2)*a^2*cos(2*x)^2 + 9*sqrt(2)*a^2*cos(x)^2 + sqrt(2)*a^2*sin(3*x)^2 + 9*sqrt(2)*a^2*sin(2*x)
^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) + 9*sqrt(2)*a^2*sin(x)^2 + 6*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2 + 2*(3*sqrt(
2)*a^2*cos(2*x) + 3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(3*x) + 6*(3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(2*
x) + 6*(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(x))*sin(3*x))*integrate(1/9*(x*cos(3/2*x)*sin(4*x) + 4*x*cos(3/
2*x)*sin(3*x) + 6*x*cos(3/2*x)*sin(2*x) - x*cos(4*x)*sin(3/2*x) - 4*x*cos(3*x)*sin(3/2*x) - 6*x*cos(2*x)*sin(3
/2*x) + 4*x*cos(3/2*x)*sin(x) - (4*x*cos(x) + x)*sin(3/2*x))/(a^2*cos(4*x)^2 + 16*a^2*cos(3*x)^2 + 36*a^2*cos(
2*x)^2 + 16*a^2*cos(x)^2 + a^2*sin(4*x)^2 + 16*a^2*sin(3*x)^2 + 36*a^2*sin(2*x)^2 + 48*a^2*sin(2*x)*sin(x) + 1
6*a^2*sin(x)^2 + 8*a^2*cos(x) + a^2 + 2*(4*a^2*cos(3*x) + 6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(4*x) + 8*(6
*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(3*x) + 12*(4*a^2*cos(x) + a^2)*cos(2*x) + 4*(2*a^2*sin(3*x) + 3*a^2*si
n(2*x) + 2*a^2*sin(x))*sin(4*x) + 16*(3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(3*x)), x) + (12*sqrt(2)*x*sin(3/2*x)
+ (9*sqrt(2)*x^2 - 8*sqrt(2))*cos(3/2*x))*sin(3*x) + 3*(12*sqrt(2)*x*sin(3/2*x) + (9*sqrt(2)*x^2 - 8*sqrt(2))*
cos(3/2*x))*sin(2*x) - (9*sqrt(2)*x^2 - 36*sqrt(2)*x*sin(x) + 3*(9*sqrt(2)*x^2 - 8*sqrt(2))*cos(x) - 8*sqrt(2)
)*sin(3/2*x))*cos(2/3*arctan2(sin(3/2*x), cos(3/2*x))) - 486*(sqrt(2)*a^2*cos(3*x)^3 + 9*sqrt(2)*a^2*cos(2*x)^
2 + 9*sqrt(2)*a^2*cos(x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) + 9*sqrt(2)*a^2*sin(x)^
2 + 6*sqrt(2)*a^2*cos(x) + 3*(2*sqrt(2)*a^2*cos(2*x) + 2*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(3*x)^2 + (sqrt(
2)*a^2*cos(3*x) + sqrt(2)*a^2)*sin(3*x)^2 + sqrt(2)*a^2 + 3*(3*sqrt(2)*a^2*cos(2*x)^2 + 3*sqrt(2)*a^2*cos(x)^2
+ 3*sqrt(2)*a^2*sin(2*x)^2 + 6*sqrt(2)*a^2*sin(2*x)*sin(x) + 3*sqrt(2)*a^2*sin(x)^2 + 4*sqrt(2)*a^2*cos(x) +
sqrt(2)*a^2 + 2*(3*sqrt(2)*a^2*cos(x) + 2*sqrt(2)*a^2)*cos(2*x))*cos(3*x) + 6*(3*sqrt(2)*a^2*cos(x) + sqrt(2)*
a^2)*cos(2*x) + 6*(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(x) + (sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(x))*cos
(3*x))*sin(3*x))*integrate(1/9*(x^2*cos(4*x)*cos(3/2*x) + 4*x^2*cos(3*x)*cos(3/2*x) + 6*x^2*cos(2*x)*cos(3/2*x
) + x^2*sin(4*x)*sin(3/2*x) + 4*x^2*sin(3*x)*sin(3/2*x) + 6*x^2*sin(2*x)*sin(3/2*x) + 4*x^2*sin(3/2*x)*sin(x)
+ (4*x^2*cos(x) + x^2)*cos(3/2*x))/(a^2*cos(4*x)^2 + 16*a^2*cos(3*x)^2 + 36*a^2*cos(2*x)^2 + 16*a^2*cos(x)^2 +
a^2*sin(4*x)^2 + 16*a^2*sin(3*x)^2 + 36*a^2*sin(2*x)^2 + 48*a^2*sin(2*x)*sin(x) + 16*a^2*sin(x)^2 + 8*a^2*cos
(x) + a^2 + 2*(4*a^2*cos(3*x) + 6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(4*x) + 8*(6*a^2*cos(2*x) + 4*a^2*cos(
x) + a^2)*cos(3*x) + 12*(4*a^2*cos(x) + a^2)*cos(2*x) + 4*(2*a^2*sin(3*x) + 3*a^2*sin(2*x) + 2*a^2*sin(x))*sin
(4*x) + 16*(3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(3*x)), x) + 648*(sqrt(2)*a^2*cos(3*x)^3 + 9*sqrt(2)*a^2*cos(2*x
)^2 + 9*sqrt(2)*a^2*cos(x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) + 9*sqrt(2)*a^2*sin(x
)^2 + 6*sqrt(2)*a^2*cos(x) + 3*(2*sqrt(2)*a^2*cos(2*x) + 2*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(3*x)^2 + (sqr
t(2)*a^2*cos(3*x) + sqrt(2)*a^2)*sin(3*x)^2 + sqrt(2)*a^2 + 3*(3*sqrt(2)*a^2*cos(2*x)^2 + 3*sqrt(2)*a^2*cos(x)
^2 + 3*sqrt(2)*a^2*sin(2*x)^2 + 6*sqrt(2)*a^2*sin(2*x)*sin(x) + 3*sqrt(2)*a^2*sin(x)^2 + 4*sqrt(2)*a^2*cos(x)
+ sqrt(2)*a^2 + 2*(3*sqrt(2)*a^2*cos(x) + 2*sqrt(2)*a^2)*cos(2*x))*cos(3*x) + 6*(3*sqrt(2)*a^2*cos(x) + sqrt(2
)*a^2)*cos(2*x) + 6*(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(x) + (sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(x))*c
os(3*x))*sin(3*x))*integrate(1/9*(x*cos(3/2*x)*sin(4*x) + 4*x*cos(3/2*x)*sin(3*x) + 6*x*cos(3/2*x)*sin(2*x) -
x*cos(4*x)*sin(3/2*x) - 4*x*cos(3*x)*sin(3/2*x) - 6*x*cos(2*x)*sin(3/2*x) + 4*x*cos(3/2*x)*sin(x) - (4*x*cos(x
) + x)*sin(3/2*x))/(a^2*cos(4*x)^2 + 16*a^2*cos(3*x)^2 + 36*a^2*cos(2*x)^2 + 16*a^2*cos(x)^2 + a^2*sin(4*x)^2
+ 16*a^2*sin(3*x)^2 + 36*a^2*sin(2*x)^2 + 48*a^2*sin(2*x)*sin(x) + 16*a^2*sin(x)^2 + 8*a^2*cos(x) + a^2 + 2*(4
*a^2*cos(3*x) + 6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(4*x) + 8*(6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(3*
x) + 12*(4*a^2*cos(x) + a^2)*cos(2*x) + 4*(2*a^2*sin(3*x) + 3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(4*x) + 16*(3*a^
2*sin(2*x) + 2*a^2*sin(x))*sin(3*x)), x) + 2*(24*sqrt(2)*sin(2*x)*sin(3/2*x) - (12*sqrt(2)*x*sin(3/2*x) + (9*s
qrt(2)*x^2 - 8*sqrt(2))*cos(3/2*x))*cos(3*x) - (9*sqrt(2)*x^2 - 8*sqrt(2))*cos(3/2*x) - 12*(sqrt(2)*x - 2*sqrt
(2)*sin(x))*sin(3/2*x))*sin(3*x) - 6*((9*sqrt(2)*x^2 - 8*sqrt(2))*cos(3/2*x) + 12*(sqrt(2)*x - 2*sqrt(2)*sin(x
))*sin(3/2*x))*sin(2*x) + 2*(27*sqrt(2)*x^2*cos(x) + 9*sqrt(2)*x^2 + 36*sqrt(2)*cos(x)^2 - 36*sqrt(2)*x*sin(x)
+ 36*sqrt(2)*sin(x)^2 - 4*sqrt(2))*sin(3/2*x) - 3*(sqrt(2)*cos(3*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(
x)^2 + sqrt(2)*sin(3*x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 2*(3*sqrt
(2)*cos(2*x) + 3*sqrt(2)*cos(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*(sqrt(2)*sin
(2*x) + sqrt(2)*sin(x))*sin(3*x) + 6*sqrt(2)*cos(x) + sqrt(2))*sin(1/3*arctan2(sin(3/2*x), cos(3/2*x))))*cos(4
/3*arctan2(sin(3/2*x), cos(3/2*x))) + 4*(72*sqrt(2)*cos(2*x)^2*sin(3/2*x) + 8*sqrt(2)*sin(3*x)^2*sin(3/2*x) +
72*sqrt(2)*sin(2*x)^2*sin(3/2*x) - 2*(12*sqrt(2)*x*cos(3/2*x) - (9*sqrt(2)*x^2 - 4*sqrt(2))*sin(3/2*x))*cos(3*
x)^2 + 2*(9*(3*sqrt(2)*x^2*sin(3/2*x) - 4*sqrt(2)*x*cos(3/2*x))*cos(2*x) - 3*(12*sqrt(2)*x*cos(x) + (9*sqrt(2)
*x^2 - 8*sqrt(2))*sin(x) + 8*sqrt(2)*x)*cos(3/2*x) - 3*(12*sqrt(2)*x*sin(3/2*x) + (9*sqrt(2)*x^2 - 8*sqrt(2))*
cos(3/2*x))*sin(2*x) + (27*sqrt(2)*x^2*cos(x) + 18*sqrt(2)*x^2 - 36*sqrt(2)*x*sin(x) - 8*sqrt(2))*sin(3/2*x))*
cos(3*x) - 18*(4*sqrt(2)*x*cos(3/2*x) - (3*sqrt(2)*x^2 + 8*sqrt(2)*cos(x))*sin(3/2*x))*cos(2*x) - 6*(12*sqrt(2
)*x*cos(x) + (9*sqrt(2)*x^2 - 8*sqrt(2))*sin(x) + 4*sqrt(2)*x)*cos(3/2*x) - 486*(sqrt(2)*a^2*cos(3*x)^3 + 9*sq
rt(2)*a^2*cos(2*x)^2 + 9*sqrt(2)*a^2*cos(x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) + 9*
sqrt(2)*a^2*sin(x)^2 + 6*sqrt(2)*a^2*cos(x) + 3*(2*sqrt(2)*a^2*cos(2*x) + 2*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*
cos(3*x)^2 + (sqrt(2)*a^2*cos(3*x) + sqrt(2)*a^2)*sin(3*x)^2 + sqrt(2)*a^2 + 3*(3*sqrt(2)*a^2*cos(2*x)^2 + 3*s
qrt(2)*a^2*cos(x)^2 + 3*sqrt(2)*a^2*sin(2*x)^2 + 6*sqrt(2)*a^2*sin(2*x)*sin(x) + 3*sqrt(2)*a^2*sin(x)^2 + 4*sq
rt(2)*a^2*cos(x) + sqrt(2)*a^2 + 2*(3*sqrt(2)*a^2*cos(x) + 2*sqrt(2)*a^2)*cos(2*x))*cos(3*x) + 6*(3*sqrt(2)*a^
2*cos(x) + sqrt(2)*a^2)*cos(2*x) + 6*(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(x) + (sqrt(2)*a^2*sin(2*x) + sqrt
(2)*a^2*sin(x))*cos(3*x))*sin(3*x))*integrate(1/9*(x^2*cos(4*x)*cos(3/2*x) + 4*x^2*cos(3*x)*cos(3/2*x) + 6*x^2
*cos(2*x)*cos(3/2*x) + x^2*sin(4*x)*sin(3/2*x) + 4*x^2*sin(3*x)*sin(3/2*x) + 6*x^2*sin(2*x)*sin(3/2*x) + 4*x^2
*sin(3/2*x)*sin(x) + (4*x^2*cos(x) + x^2)*cos(3/2*x))/(a^2*cos(4*x)^2 + 16*a^2*cos(3*x)^2 + 36*a^2*cos(2*x)^2
+ 16*a^2*cos(x)^2 + a^2*sin(4*x)^2 + 16*a^2*sin(3*x)^2 + 36*a^2*sin(2*x)^2 + 48*a^2*sin(2*x)*sin(x) + 16*a^2*s
in(x)^2 + 8*a^2*cos(x) + a^2 + 2*(4*a^2*cos(3*x) + 6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(4*x) + 8*(6*a^2*co
s(2*x) + 4*a^2*cos(x) + a^2)*cos(3*x) + 12*(4*a^2*cos(x) + a^2)*cos(2*x) + 4*(2*a^2*sin(3*x) + 3*a^2*sin(2*x)
+ 2*a^2*sin(x))*sin(4*x) + 16*(3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(3*x)), x) + 648*(sqrt(2)*a^2*cos(3*x)^3 + 9*
sqrt(2)*a^2*cos(2*x)^2 + 9*sqrt(2)*a^2*cos(x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) +
9*sqrt(2)*a^2*sin(x)^2 + 6*sqrt(2)*a^2*cos(x) + 3*(2*sqrt(2)*a^2*cos(2*x) + 2*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2
)*cos(3*x)^2 + (sqrt(2)*a^2*cos(3*x) + sqrt(2)*a^2)*sin(3*x)^2 + sqrt(2)*a^2 + 3*(3*sqrt(2)*a^2*cos(2*x)^2 + 3
*sqrt(2)*a^2*cos(x)^2 + 3*sqrt(2)*a^2*sin(2*x)^2 + 6*sqrt(2)*a^2*sin(2*x)*sin(x) + 3*sqrt(2)*a^2*sin(x)^2 + 4*
sqrt(2)*a^2*cos(x) + sqrt(2)*a^2 + 2*(3*sqrt(2)*a^2*cos(x) + 2*sqrt(2)*a^2)*cos(2*x))*cos(3*x) + 6*(3*sqrt(2)*
a^2*cos(x) + sqrt(2)*a^2)*cos(2*x) + 6*(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(x) + (sqrt(2)*a^2*sin(2*x) + sq
rt(2)*a^2*sin(x))*cos(3*x))*sin(3*x))*integrate(1/9*(x*cos(3/2*x)*sin(4*x) + 4*x*cos(3/2*x)*sin(3*x) + 6*x*cos
(3/2*x)*sin(2*x) - x*cos(4*x)*sin(3/2*x) - 4*x*cos(3*x)*sin(3/2*x) - 6*x*cos(2*x)*sin(3/2*x) + 4*x*cos(3/2*x)*
sin(x) - (4*x*cos(x) + x)*sin(3/2*x))/(a^2*cos(4*x)^2 + 16*a^2*cos(3*x)^2 + 36*a^2*cos(2*x)^2 + 16*a^2*cos(x)^
2 + a^2*sin(4*x)^2 + 16*a^2*sin(3*x)^2 + 36*a^2*sin(2*x)^2 + 48*a^2*sin(2*x)*sin(x) + 16*a^2*sin(x)^2 + 8*a^2*
cos(x) + a^2 + 2*(4*a^2*cos(3*x) + 6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(4*x) + 8*(6*a^2*cos(2*x) + 4*a^2*c
os(x) + a^2)*cos(3*x) + 12*(4*a^2*cos(x) + a^2)*cos(2*x) + 4*(2*a^2*sin(3*x) + 3*a^2*sin(2*x) + 2*a^2*sin(x))*
sin(4*x) + 16*(3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(3*x)), x) + 2*(24*sqrt(2)*sin(2*x)*sin(3/2*x) - (12*sqrt(2)*
x*sin(3/2*x) + (9*sqrt(2)*x^2 - 8*sqrt(2))*cos(3/2*x))*cos(3*x) - (9*sqrt(2)*x^2 - 8*sqrt(2))*cos(3/2*x) - 12*
(sqrt(2)*x - 2*sqrt(2)*sin(x))*sin(3/2*x))*sin(3*x) - 6*((9*sqrt(2)*x^2 - 8*sqrt(2))*cos(3/2*x) + 12*(sqrt(2)*
x - 2*sqrt(2)*sin(x))*sin(3/2*x))*sin(2*x) + 2*(27*sqrt(2)*x^2*cos(x) + 9*sqrt(2)*x^2 + 36*sqrt(2)*cos(x)^2 -
36*sqrt(2)*x*sin(x) + 36*sqrt(2)*sin(x)^2 - 4*sqrt(2))*sin(3/2*x) - 3*(sqrt(2)*cos(3*x)^2 + 9*sqrt(2)*cos(2*x)
^2 + 9*sqrt(2)*cos(x)^2 + sqrt(2)*sin(3*x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*s
in(x)^2 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*cos(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*
x) + 6*(sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*sin(3*x) + 6*sqrt(2)*cos(x) + sqrt(2))*sin(1/3*arctan2(sin(3/2*x),
cos(3/2*x))))*cos(2/3*arctan2(sin(3/2*x), cos(3/2*x))) + 4*(sqrt(2)*sin(3*x)^3 + 6*(sqrt(2)*sin(2*x) + sqrt(2)
*sin(x))*sin(3*x)^2 + (sqrt(2)*cos(3*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + 9*sqrt(2)*sin(2*x)^2 +
18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*cos(x) + sqrt(2))*cos(3*x
) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*sqrt(2)*cos(x) + sqrt(2))*sin(3*x))*cos(1/3*arctan2(sin(3/2*x)
, cos(3/2*x))) - 324*(sqrt(2)*a^2*cos(3*x)^4 + sqrt(2)*a^2*sin(3*x)^4 + 9*sqrt(2)*a^2*cos(2*x)^2 + 9*sqrt(2)*a
^2*cos(x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) + 9*sqrt(2)*a^2*sin(x)^2 + 2*(3*sqrt(2
)*a^2*cos(2*x) + 3*sqrt(2)*a^2*cos(x) + 2*sqrt(2)*a^2)*cos(3*x)^3 + 6*(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(
x))*sin(3*x)^3 + 6*sqrt(2)*a^2*cos(x) + 3*(3*sqrt(2)*a^2*cos(2*x)^2 + 3*sqrt(2)*a^2*cos(x)^2 + 3*sqrt(2)*a^2*s
in(2*x)^2 + 6*sqrt(2)*a^2*sin(2*x)*sin(x) + 3*sqrt(2)*a^2*sin(x)^2 + 6*sqrt(2)*a^2*cos(x) + 2*sqrt(2)*a^2 + 6*
(sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(2*x))*cos(3*x)^2 + (2*sqrt(2)*a^2*cos(3*x)^2 + 9*sqrt(2)*a^2*cos(2*x)^2
+ 9*sqrt(2)*a^2*cos(x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) + 9*sqrt(2)*a^2*sin(x)^2
+ 6*sqrt(2)*a^2*cos(x) + 2*sqrt(2)*a^2 + 2*(3*sqrt(2)*a^2*cos(2*x) + 3*sqrt(2)*a^2*cos(x) + 2*sqrt(2)*a^2)*co
s(3*x) + 6*(3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(2*x))*sin(3*x)^2 + sqrt(2)*a^2 + 2*(9*sqrt(2)*a^2*cos(2*x)
^2 + 9*sqrt(2)*a^2*cos(x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) + 9*sqrt(2)*a^2*sin(x)
^2 + 9*sqrt(2)*a^2*cos(x) + 2*sqrt(2)*a^2 + 9*(2*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(2*x))*cos(3*x) + 6*(3*s
qrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(2*x) + 6*(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(x) + (sqrt(2)*a^2*sin(2*
x) + sqrt(2)*a^2*sin(x))*cos(3*x)^2 + 2*(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(x))*cos(3*x))*sin(3*x))*integr
ate(1/9*(x^2*cos(4*x)*cos(3/2*x) + 4*x^2*cos(3*x)*cos(3/2*x) + 6*x^2*cos(2*x)*cos(3/2*x) + x^2*sin(4*x)*sin(3/
2*x) + 4*x^2*sin(3*x)*sin(3/2*x) + 6*x^2*sin(2*x)*sin(3/2*x) + 4*x^2*sin(3/2*x)*sin(x) + (4*x^2*cos(x) + x^2)*
cos(3/2*x))/(a^2*cos(4*x)^2 + 16*a^2*cos(3*x)^2 + 36*a^2*cos(2*x)^2 + 16*a^2*cos(x)^2 + a^2*sin(4*x)^2 + 16*a^
2*sin(3*x)^2 + 36*a^2*sin(2*x)^2 + 48*a^2*sin(2*x)*sin(x) + 16*a^2*sin(x)^2 + 8*a^2*cos(x) + a^2 + 2*(4*a^2*co
s(3*x) + 6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(4*x) + 8*(6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(3*x) + 12
*(4*a^2*cos(x) + a^2)*cos(2*x) + 4*(2*a^2*sin(3*x) + 3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(4*x) + 16*(3*a^2*sin(2
*x) + 2*a^2*sin(x))*sin(3*x)), x) + 432*(sqrt(2)*a^2*cos(3*x)^4 + sqrt(2)*a^2*sin(3*x)^4 + 9*sqrt(2)*a^2*cos(2
*x)^2 + 9*sqrt(2)*a^2*cos(x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) + 9*sqrt(2)*a^2*sin
(x)^2 + 2*(3*sqrt(2)*a^2*cos(2*x) + 3*sqrt(2)*a^2*cos(x) + 2*sqrt(2)*a^2)*cos(3*x)^3 + 6*(sqrt(2)*a^2*sin(2*x)
+ sqrt(2)*a^2*sin(x))*sin(3*x)^3 + 6*sqrt(2)*a^2*cos(x) + 3*(3*sqrt(2)*a^2*cos(2*x)^2 + 3*sqrt(2)*a^2*cos(x)^
2 + 3*sqrt(2)*a^2*sin(2*x)^2 + 6*sqrt(2)*a^2*sin(2*x)*sin(x) + 3*sqrt(2)*a^2*sin(x)^2 + 6*sqrt(2)*a^2*cos(x) +
2*sqrt(2)*a^2 + 6*(sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(2*x))*cos(3*x)^2 + (2*sqrt(2)*a^2*cos(3*x)^2 + 9*sqr
t(2)*a^2*cos(2*x)^2 + 9*sqrt(2)*a^2*cos(x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) + 9*s
qrt(2)*a^2*sin(x)^2 + 6*sqrt(2)*a^2*cos(x) + 2*sqrt(2)*a^2 + 2*(3*sqrt(2)*a^2*cos(2*x) + 3*sqrt(2)*a^2*cos(x)
+ 2*sqrt(2)*a^2)*cos(3*x) + 6*(3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(2*x))*sin(3*x)^2 + sqrt(2)*a^2 + 2*(9*s
qrt(2)*a^2*cos(2*x)^2 + 9*sqrt(2)*a^2*cos(x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) + 9
*sqrt(2)*a^2*sin(x)^2 + 9*sqrt(2)*a^2*cos(x) + 2*sqrt(2)*a^2 + 9*(2*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(2*x)
)*cos(3*x) + 6*(3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(2*x) + 6*(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(x) +
(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(x))*cos(3*x)^2 + 2*(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(x))*cos(3*x
))*sin(3*x))*integrate(1/9*(x*cos(3/2*x)*sin(4*x) + 4*x*cos(3/2*x)*sin(3*x) + 6*x*cos(3/2*x)*sin(2*x) - x*cos(
4*x)*sin(3/2*x) - 4*x*cos(3*x)*sin(3/2*x) - 6*x*cos(2*x)*sin(3/2*x) + 4*x*cos(3/2*x)*sin(x) - (4*x*cos(x) + x)
*sin(3/2*x))/(a^2*cos(4*x)^2 + 16*a^2*cos(3*x)^2 + 36*a^2*cos(2*x)^2 + 16*a^2*cos(x)^2 + a^2*sin(4*x)^2 + 16*a
^2*sin(3*x)^2 + 36*a^2*sin(2*x)^2 + 48*a^2*sin(2*x)*sin(x) + 16*a^2*sin(x)^2 + 8*a^2*cos(x) + a^2 + 2*(4*a^2*c
os(3*x) + 6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(4*x) + 8*(6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(3*x) + 1
2*(4*a^2*cos(x) + a^2)*cos(2*x) + 4*(2*a^2*sin(3*x) + 3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(4*x) + 16*(3*a^2*sin(
2*x) + 2*a^2*sin(x))*sin(3*x)), x) + (sqrt(2)*cos(3*x)^4 + sqrt(2)*sin(3*x)^4 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt
(2)*cos(x) + 2*sqrt(2))*cos(3*x)^3 + 6*(sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*sin(3*x)^3 + 3*(3*sqrt(2)*cos(2*x)^
2 + 3*sqrt(2)*cos(x)^2 + 3*sqrt(2)*sin(2*x)^2 + 6*sqrt(2)*sin(2*x)*sin(x) + 3*sqrt(2)*sin(x)^2 + 6*(sqrt(2)*co
s(x) + sqrt(2))*cos(2*x) + 6*sqrt(2)*cos(x) + 2*sqrt(2))*cos(3*x)^2 + 9*(sqrt(2)*cos(3*x)^2 + 9*sqrt(2)*cos(2*
x)^2 + 9*sqrt(2)*cos(x)^2 + sqrt(2)*sin(3*x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)
*sin(x)^2 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*cos(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(
2*x) + 6*(sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*sin(3*x) + 6*sqrt(2)*cos(x) + sqrt(2))*cos(4/3*arctan2(sin(3/2*x)
, cos(3/2*x)))^2 + 9*(sqrt(2)*cos(3*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + sqrt(2)*sin(3*x)^2 + 9*
sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*cos(x
) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*(sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*sin(3*
x) + 6*sqrt(2)*cos(x) + sqrt(2))*cos(2/3*arctan2(sin(3/2*x), cos(3/2*x)))^2 + (2*sqrt(2)*cos(3*x)^2 + 9*sqrt(2
)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 2
*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*cos(x) + 2*sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*sq
rt(2)*cos(x) + 2*sqrt(2))*sin(3*x)^2 + 9*(sqrt(2)*cos(3*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + sqr
t(2)*sin(3*x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(2*
x) + 3*sqrt(2)*cos(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*(sqrt(2)*sin(2*x) + sq
rt(2)*sin(x))*sin(3*x) + 6*sqrt(2)*cos(x) + sqrt(2))*sin(4/3*arctan2(sin(3/2*x), cos(3/2*x)))^2 + 9*(sqrt(2)*c
os(3*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + sqrt(2)*sin(3*x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)
*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*cos(x) + sqrt(2))*cos(3*x) + 6*(3*sq
rt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*(sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*sin(3*x) + 6*sqrt(2)*cos(x) + sqrt(2)
)*sin(2/3*arctan2(sin(3/2*x), cos(3/2*x)))^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + 9*sqrt(2)*sin(2*x)^
2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 2*(9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + 9*sqrt(2)
*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 9*(2*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 9*sq
rt(2)*cos(x) + 2*sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*(sqrt(2)*cos(3*x)^3 + 3*(2*sq
rt(2)*cos(2*x) + 2*sqrt(2)*cos(x) + sqrt(2))*cos(3*x)^2 + (sqrt(2)*cos(3*x) + sqrt(2))*sin(3*x)^2 + 9*sqrt(2)*
cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 3*(
3*sqrt(2)*cos(2*x)^2 + 3*sqrt(2)*cos(x)^2 + 3*sqrt(2)*sin(2*x)^2 + 6*sqrt(2)*sin(2*x)*sin(x) + 3*sqrt(2)*sin(x
)^2 + 2*(3*sqrt(2)*cos(x) + 2*sqrt(2))*cos(2*x) + 4*sqrt(2)*cos(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) +
sqrt(2))*cos(2*x) + 3*(sqrt(2)*cos(3*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + sqrt(2)*sin(3*x)^2 +
9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*cos
(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*(sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*sin(
3*x) + 6*sqrt(2)*cos(x) + sqrt(2))*cos(2/3*arctan2(sin(3/2*x), cos(3/2*x))) + 6*((sqrt(2)*sin(2*x) + sqrt(2)*s
in(x))*cos(3*x) + sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*sin(3*x) + 6*sqrt(2)*cos(x) + sqrt(2))*cos(4/3*arctan2(si
n(3/2*x), cos(3/2*x))) + 6*(sqrt(2)*cos(3*x)^3 + 3*(2*sqrt(2)*cos(2*x) + 2*sqrt(2)*cos(x) + sqrt(2))*cos(3*x)^
2 + (sqrt(2)*cos(3*x) + sqrt(2))*sin(3*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + 9*sqrt(2)*sin(2*x)^2
+ 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 3*(3*sqrt(2)*cos(2*x)^2 + 3*sqrt(2)*cos(x)^2 + 3*sqrt(2)*
sin(2*x)^2 + 6*sqrt(2)*sin(2*x)*sin(x) + 3*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(x) + 2*sqrt(2))*cos(2*x) + 4*sq
rt(2)*cos(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*((sqrt(2)*sin(2*x) + sqrt(2)*si
n(x))*cos(3*x) + sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*sin(3*x) + 6*sqrt(2)*cos(x) + sqrt(2))*cos(2/3*arctan2(sin
(3/2*x), cos(3/2*x))) + 6*((sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*cos(3*x)^2 + 2*(sqrt(2)*sin(2*x) + sqrt(2)*sin(
x))*cos(3*x) + sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*sin(3*x) + 6*(sqrt(2)*sin(3*x)^3 + 6*(sqrt(2)*sin(2*x) + sqr
t(2)*sin(x))*sin(3*x)^2 + (sqrt(2)*cos(3*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + 9*sqrt(2)*sin(2*x)
^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*cos(x) + sqrt(2))*cos
(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*sqrt(2)*cos(x) + sqrt(2))*sin(3*x) + 3*(sqrt(2)*cos(3*x)^2
+ 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + sqrt(2)*sin(3*x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)
*sin(x) + 9*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*cos(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos
(x) + sqrt(2))*cos(2*x) + 6*(sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*sin(3*x) + 6*sqrt(2)*cos(x) + sqrt(2))*sin(2/3
*arctan2(sin(3/2*x), cos(3/2*x))))*sin(4/3*arctan2(sin(3/2*x), cos(3/2*x))) + 6*(sqrt(2)*sin(3*x)^3 + 6*(sqrt(
2)*sin(2*x) + sqrt(2)*sin(x))*sin(3*x)^2 + (sqrt(2)*cos(3*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + 9
*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*cos(
x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*sqrt(2)*cos(x) + sqrt(2))*sin(3*x))*sin(2
/3*arctan2(sin(3/2*x), cos(3/2*x))) + 6*sqrt(2)*cos(x) + sqrt(2))*log(cos(1/3*arctan2(sin(3/2*x), cos(3/2*x)))
^2 + sin(1/3*arctan2(sin(3/2*x), cos(3/2*x)))^2 + 2*sin(1/3*arctan2(sin(3/2*x), cos(3/2*x))) + 1) - (sqrt(2)*c
os(3*x)^4 + sqrt(2)*sin(3*x)^4 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*cos(x) + 2*sqrt(2))*cos(3*x)^3 + 6*(sqrt(2)
*sin(2*x) + sqrt(2)*sin(x))*sin(3*x)^3 + 3*(3*sqrt(2)*cos(2*x)^2 + 3*sqrt(2)*cos(x)^2 + 3*sqrt(2)*sin(2*x)^2 +
6*sqrt(2)*sin(2*x)*sin(x) + 3*sqrt(2)*sin(x)^2 + 6*(sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*sqrt(2)*cos(x) + 2
*sqrt(2))*cos(3*x)^2 + 9*(sqrt(2)*cos(3*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + sqrt(2)*sin(3*x)^2
+ 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*c
os(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*(sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*si
n(3*x) + 6*sqrt(2)*cos(x) + sqrt(2))*cos(4/3*arctan2(sin(3/2*x), cos(3/2*x)))^2 + 9*(sqrt(2)*cos(3*x)^2 + 9*sq
rt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + sqrt(2)*sin(3*x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x)
+ 9*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*cos(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) + s
qrt(2))*cos(2*x) + 6*(sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*sin(3*x) + 6*sqrt(2)*cos(x) + sqrt(2))*cos(2/3*arctan
2(sin(3/2*x), cos(3/2*x)))^2 + (2*sqrt(2)*cos(3*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + 9*sqrt(2)*s
in(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*cos(x) + 2*sqr
t(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*sqrt(2)*cos(x) + 2*sqrt(2))*sin(3*x)^2 + 9*(sqrt(
2)*cos(3*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + sqrt(2)*sin(3*x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqr
t(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*cos(x) + sqrt(2))*cos(3*x) + 6*(
3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*(sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*sin(3*x) + 6*sqrt(2)*cos(x) + sqr
t(2))*sin(4/3*arctan2(sin(3/2*x), cos(3/2*x)))^2 + 9*(sqrt(2)*cos(3*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*co
s(x)^2 + sqrt(2)*sin(3*x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 2*(3*sq
rt(2)*cos(2*x) + 3*sqrt(2)*cos(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*(sqrt(2)*s
in(2*x) + sqrt(2)*sin(x))*sin(3*x) + 6*sqrt(2)*cos(x) + sqrt(2))*sin(2/3*arctan2(sin(3/2*x), cos(3/2*x)))^2 +
9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(
x)^2 + 2*(9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sq
rt(2)*sin(x)^2 + 9*(2*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 9*sqrt(2)*cos(x) + 2*sqrt(2))*cos(3*x) + 6*(3*sqrt(
2)*cos(x) + sqrt(2))*cos(2*x) + 6*(sqrt(2)*cos(3*x)^3 + 3*(2*sqrt(2)*cos(2*x) + 2*sqrt(2)*cos(x) + sqrt(2))*co
s(3*x)^2 + (sqrt(2)*cos(3*x) + sqrt(2))*sin(3*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + 9*sqrt(2)*sin
(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 3*(3*sqrt(2)*cos(2*x)^2 + 3*sqrt(2)*cos(x)^2 + 3*s
qrt(2)*sin(2*x)^2 + 6*sqrt(2)*sin(2*x)*sin(x) + 3*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(x) + 2*sqrt(2))*cos(2*x)
+ 4*sqrt(2)*cos(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 3*(sqrt(2)*cos(3*x)^2 + 9*
sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + sqrt(2)*sin(3*x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(
x) + 9*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*cos(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) +
sqrt(2))*cos(2*x) + 6*(sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*sin(3*x) + 6*sqrt(2)*cos(x) + sqrt(2))*cos(2/3*arct
an2(sin(3/2*x), cos(3/2*x))) + 6*((sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*cos(3*x) + sqrt(2)*sin(2*x) + sqrt(2)*si
n(x))*sin(3*x) + 6*sqrt(2)*cos(x) + sqrt(2))*cos(4/3*arctan2(sin(3/2*x), cos(3/2*x))) + 6*(sqrt(2)*cos(3*x)^3
+ 3*(2*sqrt(2)*cos(2*x) + 2*sqrt(2)*cos(x) + sqrt(2))*cos(3*x)^2 + (sqrt(2)*cos(3*x) + sqrt(2))*sin(3*x)^2 + 9
*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x
)^2 + 3*(3*sqrt(2)*cos(2*x)^2 + 3*sqrt(2)*cos(x)^2 + 3*sqrt(2)*sin(2*x)^2 + 6*sqrt(2)*sin(2*x)*sin(x) + 3*sqrt
(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(x) + 2*sqrt(2))*cos(2*x) + 4*sqrt(2)*cos(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)
*cos(x) + sqrt(2))*cos(2*x) + 6*((sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*cos(3*x) + sqrt(2)*sin(2*x) + sqrt(2)*sin
(x))*sin(3*x) + 6*sqrt(2)*cos(x) + sqrt(2))*cos(2/3*arctan2(sin(3/2*x), cos(3/2*x))) + 6*((sqrt(2)*sin(2*x) +
sqrt(2)*sin(x))*cos(3*x)^2 + 2*(sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*cos(3*x) + sqrt(2)*sin(2*x) + sqrt(2)*sin(x
))*sin(3*x) + 6*(sqrt(2)*sin(3*x)^3 + 6*(sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*sin(3*x)^2 + (sqrt(2)*cos(3*x)^2 +
9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin
(x)^2 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*cos(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x)
+ 6*sqrt(2)*cos(x) + sqrt(2))*sin(3*x) + 3*(sqrt(2)*cos(3*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 +
sqrt(2)*sin(3*x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos
(2*x) + 3*sqrt(2)*cos(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*(sqrt(2)*sin(2*x) +
sqrt(2)*sin(x))*sin(3*x) + 6*sqrt(2)*cos(x) + sqrt(2))*sin(2/3*arctan2(sin(3/2*x), cos(3/2*x))))*sin(4/3*arct
an2(sin(3/2*x), cos(3/2*x))) + 6*(sqrt(2)*sin(3*x)^3 + 6*(sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*sin(3*x)^2 + (sqr
t(2)*cos(3*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x
) + 9*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*cos(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) +
sqrt(2))*cos(2*x) + 6*sqrt(2)*cos(x) + sqrt(2))*sin(3*x))*sin(2/3*arctan2(sin(3/2*x), cos(3/2*x))) + 6*sqrt(2)
*cos(x) + sqrt(2))*log(cos(1/3*arctan2(sin(3/2*x), cos(3/2*x)))^2 + sin(1/3*arctan2(sin(3/2*x), cos(3/2*x)))^2
- 2*sin(1/3*arctan2(sin(3/2*x), cos(3/2*x))) + 1) - 4*(24*sqrt(2)*cos(2*x)^2*cos(3/2*x) + 24*sqrt(2)*cos(3/2*
x)*sin(2*x)^2 + (3*sqrt(2)*x^2*cos(3/2*x) + 4*sqrt(2)*x*sin(3/2*x))*cos(3*x)^2 + 16*(3*sqrt(2)*cos(x) + sqrt(2
))*cos(2*x)*cos(3/2*x) + 2*(8*sqrt(2)*cos(2*x)*cos(3/2*x) - 8*sqrt(2)*sin(2*x)*sin(3/2*x) + (3*sqrt(2)*x^2 + 8
*sqrt(2)*cos(x))*cos(3/2*x) + 4*(sqrt(2)*x - 2*sqrt(2)*sin(x))*sin(3/2*x))*cos(3*x) + (3*sqrt(2)*x^2 + 24*sqrt
(2)*cos(x)^2 + 24*sqrt(2)*sin(x)^2 + 16*sqrt(2)*cos(x))*cos(3/2*x) + 16*(3*sqrt(2)*cos(3/2*x)*sin(x) - sqrt(2)
*sin(3/2*x))*sin(2*x) + 4*(sqrt(2)*x - 4*sqrt(2)*sin(x))*sin(3/2*x))*sin(3*x) - 4*((9*sqrt(2)*x^2 - 8*sqrt(2))
*cos(3/2*x) + 12*(sqrt(2)*x - 4*sqrt(2)*sin(x))*sin(3/2*x))*sin(2*x) + 4*(3*sqrt(2)*x^2 + 24*sqrt(2)*cos(x)^2
- 12*sqrt(2)*x*sin(x) + 24*sqrt(2)*sin(x)^2 + (9*sqrt(2)*x^2 + 8*sqrt(2))*cos(x))*sin(3/2*x) + 4*(sqrt(2)*cos(
3*x)^3 + 3*(2*sqrt(2)*cos(2*x) + 2*sqrt(2)*cos(x) + sqrt(2))*cos(3*x)^2 + (sqrt(2)*cos(3*x) + sqrt(2))*sin(3*x
)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2
)*sin(x)^2 + 3*(3*sqrt(2)*cos(2*x)^2 + 3*sqrt(2)*cos(x)^2 + 3*sqrt(2)*sin(2*x)^2 + 6*sqrt(2)*sin(2*x)*sin(x) +
3*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(x) + 2*sqrt(2))*cos(2*x) + 4*sqrt(2)*cos(x) + sqrt(2))*cos(3*x) + 6*(3*
sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 3*(sqrt(2)*cos(3*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + sqrt(
2)*sin(3*x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(2*x)
+ 3*sqrt(2)*cos(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*(sqrt(2)*sin(2*x) + sqrt
(2)*sin(x))*sin(3*x) + 6*sqrt(2)*cos(x) + sqrt(2))*cos(4/3*arctan2(sin(3/2*x), cos(3/2*x))) + 3*(sqrt(2)*cos(3
*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + sqrt(2)*sin(3*x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin
(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*cos(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2
)*cos(x) + sqrt(2))*cos(2*x) + 6*(sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*sin(3*x) + 6*sqrt(2)*cos(x) + sqrt(2))*co
s(2/3*arctan2(sin(3/2*x), cos(3/2*x))) + 6*((sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*cos(3*x) + sqrt(2)*sin(2*x) +
sqrt(2)*sin(x))*sin(3*x) + 6*sqrt(2)*cos(x) + sqrt(2))*sin(5/3*arctan2(sin(3/2*x), cos(3/2*x))) - 4*(8*sqrt(2)
*cos(3*x)^2*cos(3/2*x) + 72*sqrt(2)*cos(2*x)^2*cos(3/2*x) + 72*sqrt(2)*cos(3/2*x)*sin(2*x)^2 + 144*sqrt(2)*cos
(3/2*x)*sin(2*x)*sin(x) + 48*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x)*cos(3/2*x) + 2*(12*sqrt(2)*x*sin(3/2*x) + (
9*sqrt(2)*x^2 - 4*sqrt(2))*cos(3/2*x))*sin(3*x)^2 + 16*(3*sqrt(2)*cos(2*x)*cos(3/2*x) + (3*sqrt(2)*cos(x) + sq
rt(2))*cos(3/2*x))*cos(3*x) + 8*(9*sqrt(2)*cos(x)^2 + 9*sqrt(2)*sin(x)^2 + 6*sqrt(2)*cos(x) + sqrt(2))*cos(3/2
*x) - 3*(sqrt(2)*cos(3*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + sqrt(2)*sin(3*x)^2 + 9*sqrt(2)*sin(2
*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*cos(x) + sqrt(2))*
cos(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*(sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*sin(3*x) + 6*sqrt(2
)*cos(x) + sqrt(2))*cos(1/3*arctan2(sin(3/2*x), cos(3/2*x))) + 486*(sqrt(2)*a^2*sin(3*x)^3 + 6*(sqrt(2)*a^2*si
n(2*x) + sqrt(2)*a^2*sin(x))*sin(3*x)^2 + (sqrt(2)*a^2*cos(3*x)^2 + 9*sqrt(2)*a^2*cos(2*x)^2 + 9*sqrt(2)*a^2*c
os(x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) + 9*sqrt(2)*a^2*sin(x)^2 + 6*sqrt(2)*a^2*c
os(x) + sqrt(2)*a^2 + 2*(3*sqrt(2)*a^2*cos(2*x) + 3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(3*x) + 6*(3*sqrt(2)*
a^2*cos(x) + sqrt(2)*a^2)*cos(2*x))*sin(3*x))*integrate(1/9*(x^2*cos(4*x)*cos(3/2*x) + 4*x^2*cos(3*x)*cos(3/2*
x) + 6*x^2*cos(2*x)*cos(3/2*x) + x^2*sin(4*x)*sin(3/2*x) + 4*x^2*sin(3*x)*sin(3/2*x) + 6*x^2*sin(2*x)*sin(3/2*
x) + 4*x^2*sin(3/2*x)*sin(x) + (4*x^2*cos(x) + x^2)*cos(3/2*x))/(a^2*cos(4*x)^2 + 16*a^2*cos(3*x)^2 + 36*a^2*c
os(2*x)^2 + 16*a^2*cos(x)^2 + a^2*sin(4*x)^2 + 16*a^2*sin(3*x)^2 + 36*a^2*sin(2*x)^2 + 48*a^2*sin(2*x)*sin(x)
+ 16*a^2*sin(x)^2 + 8*a^2*cos(x) + a^2 + 2*(4*a^2*cos(3*x) + 6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(4*x) + 8
*(6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(3*x) + 12*(4*a^2*cos(x) + a^2)*cos(2*x) + 4*(2*a^2*sin(3*x) + 3*a^2
*sin(2*x) + 2*a^2*sin(x))*sin(4*x) + 16*(3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(3*x)), x) - 648*(sqrt(2)*a^2*sin(3
*x)^3 + 6*(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(x))*sin(3*x)^2 + (sqrt(2)*a^2*cos(3*x)^2 + 9*sqrt(2)*a^2*cos
(2*x)^2 + 9*sqrt(2)*a^2*cos(x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) + 9*sqrt(2)*a^2*s
in(x)^2 + 6*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2 + 2*(3*sqrt(2)*a^2*cos(2*x) + 3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)
*cos(3*x) + 6*(3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(2*x))*sin(3*x))*integrate(1/9*(x*cos(3/2*x)*sin(4*x) +
4*x*cos(3/2*x)*sin(3*x) + 6*x*cos(3/2*x)*sin(2*x) - x*cos(4*x)*sin(3/2*x) - 4*x*cos(3*x)*sin(3/2*x) - 6*x*cos(
2*x)*sin(3/2*x) + 4*x*cos(3/2*x)*sin(x) - (4*x*cos(x) + x)*sin(3/2*x))/(a^2*cos(4*x)^2 + 16*a^2*cos(3*x)^2 + 3
6*a^2*cos(2*x)^2 + 16*a^2*cos(x)^2 + a^2*sin(4*x)^2 + 16*a^2*sin(3*x)^2 + 36*a^2*sin(2*x)^2 + 48*a^2*sin(2*x)*
sin(x) + 16*a^2*sin(x)^2 + 8*a^2*cos(x) + a^2 + 2*(4*a^2*cos(3*x) + 6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(4
*x) + 8*(6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(3*x) + 12*(4*a^2*cos(x) + a^2)*cos(2*x) + 4*(2*a^2*sin(3*x)
+ 3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(4*x) + 16*(3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(3*x)), x) + 2*((12*sqrt(2)*
x*cos(3/2*x) - (9*sqrt(2)*x^2 - 8*sqrt(2))*sin(3/2*x))*cos(3*x) + 3*(12*sqrt(2)*x*cos(3/2*x) - (9*sqrt(2)*x^2
- 8*sqrt(2))*sin(3/2*x))*cos(2*x) + 3*(9*sqrt(2)*x^2*sin(x) + 12*sqrt(2)*x*cos(x) + 4*sqrt(2)*x)*cos(3/2*x) +
9*(3*sqrt(2)*x^2*cos(3/2*x) + 4*sqrt(2)*x*sin(3/2*x))*sin(2*x) - (9*sqrt(2)*x^2 - 36*sqrt(2)*x*sin(x) + 3*(9*s
qrt(2)*x^2 - 8*sqrt(2))*cos(x) - 8*sqrt(2))*sin(3/2*x))*sin(3*x) + 6*((12*sqrt(2)*x*cos(3/2*x) - (9*sqrt(2)*x^
2 - 8*sqrt(2))*sin(3/2*x))*cos(3*x) + 3*(12*sqrt(2)*x*cos(3/2*x) - (9*sqrt(2)*x^2 - 8*sqrt(2))*sin(3/2*x))*cos
(2*x) + 3*(12*sqrt(2)*x*cos(x) + (9*sqrt(2)*x^2 - 8*sqrt(2))*sin(x) + 4*sqrt(2)*x)*cos(3/2*x) + 243*(sqrt(2)*a
^2*cos(3*x)^2 + 9*sqrt(2)*a^2*cos(2*x)^2 + 9*sqrt(2)*a^2*cos(x)^2 + sqrt(2)*a^2*sin(3*x)^2 + 9*sqrt(2)*a^2*sin
(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) + 9*sqrt(2)*a^2*sin(x)^2 + 6*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2 + 2*(3*
sqrt(2)*a^2*cos(2*x) + 3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(3*x) + 6*(3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*c
os(2*x) + 6*(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(x))*sin(3*x))*integrate(1/9*(x^2*cos(4*x)*cos(3/2*x) + 4*x
^2*cos(3*x)*cos(3/2*x) + 6*x^2*cos(2*x)*cos(3/2*x) + x^2*sin(4*x)*sin(3/2*x) + 4*x^2*sin(3*x)*sin(3/2*x) + 6*x
^2*sin(2*x)*sin(3/2*x) + 4*x^2*sin(3/2*x)*sin(x) + (4*x^2*cos(x) + x^2)*cos(3/2*x))/(a^2*cos(4*x)^2 + 16*a^2*c
os(3*x)^2 + 36*a^2*cos(2*x)^2 + 16*a^2*cos(x)^2 + a^2*sin(4*x)^2 + 16*a^2*sin(3*x)^2 + 36*a^2*sin(2*x)^2 + 48*
a^2*sin(2*x)*sin(x) + 16*a^2*sin(x)^2 + 8*a^2*cos(x) + a^2 + 2*(4*a^2*cos(3*x) + 6*a^2*cos(2*x) + 4*a^2*cos(x)
+ a^2)*cos(4*x) + 8*(6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(3*x) + 12*(4*a^2*cos(x) + a^2)*cos(2*x) + 4*(2*
a^2*sin(3*x) + 3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(4*x) + 16*(3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(3*x)), x) - 32
4*(sqrt(2)*a^2*cos(3*x)^2 + 9*sqrt(2)*a^2*cos(2*x)^2 + 9*sqrt(2)*a^2*cos(x)^2 + sqrt(2)*a^2*sin(3*x)^2 + 9*sqr
t(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) + 9*sqrt(2)*a^2*sin(x)^2 + 6*sqrt(2)*a^2*cos(x) + sqrt(2)
*a^2 + 2*(3*sqrt(2)*a^2*cos(2*x) + 3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(3*x) + 6*(3*sqrt(2)*a^2*cos(x) + sq
rt(2)*a^2)*cos(2*x) + 6*(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(x))*sin(3*x))*integrate(1/9*(x*cos(3/2*x)*sin(
4*x) + 4*x*cos(3/2*x)*sin(3*x) + 6*x*cos(3/2*x)*sin(2*x) - x*cos(4*x)*sin(3/2*x) - 4*x*cos(3*x)*sin(3/2*x) - 6
*x*cos(2*x)*sin(3/2*x) + 4*x*cos(3/2*x)*sin(x) - (4*x*cos(x) + x)*sin(3/2*x))/(a^2*cos(4*x)^2 + 16*a^2*cos(3*x
)^2 + 36*a^2*cos(2*x)^2 + 16*a^2*cos(x)^2 + a^2*sin(4*x)^2 + 16*a^2*sin(3*x)^2 + 36*a^2*sin(2*x)^2 + 48*a^2*si
n(2*x)*sin(x) + 16*a^2*sin(x)^2 + 8*a^2*cos(x) + a^2 + 2*(4*a^2*cos(3*x) + 6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2
)*cos(4*x) + 8*(6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(3*x) + 12*(4*a^2*cos(x) + a^2)*cos(2*x) + 4*(2*a^2*si
n(3*x) + 3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(4*x) + 16*(3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(3*x)), x) + (12*sqrt
(2)*x*sin(3/2*x) + (9*sqrt(2)*x^2 - 8*sqrt(2))*cos(3/2*x))*sin(3*x) + 3*(12*sqrt(2)*x*sin(3/2*x) + (9*sqrt(2)*
x^2 - 8*sqrt(2))*cos(3/2*x))*sin(2*x) - (9*sqrt(2)*x^2 - 36*sqrt(2)*x*sin(x) + 3*(9*sqrt(2)*x^2 - 8*sqrt(2))*c
os(x) - 8*sqrt(2))*sin(3/2*x))*sin(2/3*arctan2(sin(3/2*x), cos(3/2*x))))*sin(4/3*arctan2(sin(3/2*x), cos(3/2*x
))) - 4*(8*sqrt(2)*cos(3*x)^2*cos(3/2*x) + 72*sqrt(2)*cos(2*x)^2*cos(3/2*x) + 72*sqrt(2)*cos(3/2*x)*sin(2*x)^2
+ 144*sqrt(2)*cos(3/2*x)*sin(2*x)*sin(x) + 48*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x)*cos(3/2*x) + 2*(12*sqrt(2
)*x*sin(3/2*x) + (9*sqrt(2)*x^2 - 4*sqrt(2))*cos(3/2*x))*sin(3*x)^2 + 16*(3*sqrt(2)*cos(2*x)*cos(3/2*x) + (3*s
qrt(2)*cos(x) + sqrt(2))*cos(3/2*x))*cos(3*x) + 8*(9*sqrt(2)*cos(x)^2 + 9*sqrt(2)*sin(x)^2 + 6*sqrt(2)*cos(x)
+ sqrt(2))*cos(3/2*x) - 3*(sqrt(2)*cos(3*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + sqrt(2)*sin(3*x)^2
+ 9*sqrt(2)*sin(2*x)^2 + 18*sqrt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(2*x) + 3*sqrt(2)*
cos(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*(sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*s
in(3*x) + 6*sqrt(2)*cos(x) + sqrt(2))*cos(1/3*arctan2(sin(3/2*x), cos(3/2*x))) + 486*(sqrt(2)*a^2*sin(3*x)^3 +
6*(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(x))*sin(3*x)^2 + (sqrt(2)*a^2*cos(3*x)^2 + 9*sqrt(2)*a^2*cos(2*x)^2
+ 9*sqrt(2)*a^2*cos(x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x) + 9*sqrt(2)*a^2*sin(x)^2
+ 6*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2 + 2*(3*sqrt(2)*a^2*cos(2*x) + 3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(3*
x) + 6*(3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(2*x))*sin(3*x))*integrate(1/9*(x^2*cos(4*x)*cos(3/2*x) + 4*x^2
*cos(3*x)*cos(3/2*x) + 6*x^2*cos(2*x)*cos(3/2*x) + x^2*sin(4*x)*sin(3/2*x) + 4*x^2*sin(3*x)*sin(3/2*x) + 6*x^2
*sin(2*x)*sin(3/2*x) + 4*x^2*sin(3/2*x)*sin(x) + (4*x^2*cos(x) + x^2)*cos(3/2*x))/(a^2*cos(4*x)^2 + 16*a^2*cos
(3*x)^2 + 36*a^2*cos(2*x)^2 + 16*a^2*cos(x)^2 + a^2*sin(4*x)^2 + 16*a^2*sin(3*x)^2 + 36*a^2*sin(2*x)^2 + 48*a^
2*sin(2*x)*sin(x) + 16*a^2*sin(x)^2 + 8*a^2*cos(x) + a^2 + 2*(4*a^2*cos(3*x) + 6*a^2*cos(2*x) + 4*a^2*cos(x) +
a^2)*cos(4*x) + 8*(6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(3*x) + 12*(4*a^2*cos(x) + a^2)*cos(2*x) + 4*(2*a^
2*sin(3*x) + 3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(4*x) + 16*(3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(3*x)), x) - 648*
(sqrt(2)*a^2*sin(3*x)^3 + 6*(sqrt(2)*a^2*sin(2*x) + sqrt(2)*a^2*sin(x))*sin(3*x)^2 + (sqrt(2)*a^2*cos(3*x)^2 +
9*sqrt(2)*a^2*cos(2*x)^2 + 9*sqrt(2)*a^2*cos(x)^2 + 9*sqrt(2)*a^2*sin(2*x)^2 + 18*sqrt(2)*a^2*sin(2*x)*sin(x)
+ 9*sqrt(2)*a^2*sin(x)^2 + 6*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2 + 2*(3*sqrt(2)*a^2*cos(2*x) + 3*sqrt(2)*a^2*cos
(x) + sqrt(2)*a^2)*cos(3*x) + 6*(3*sqrt(2)*a^2*cos(x) + sqrt(2)*a^2)*cos(2*x))*sin(3*x))*integrate(1/9*(x*cos(
3/2*x)*sin(4*x) + 4*x*cos(3/2*x)*sin(3*x) + 6*x*cos(3/2*x)*sin(2*x) - x*cos(4*x)*sin(3/2*x) - 4*x*cos(3*x)*sin
(3/2*x) - 6*x*cos(2*x)*sin(3/2*x) + 4*x*cos(3/2*x)*sin(x) - (4*x*cos(x) + x)*sin(3/2*x))/(a^2*cos(4*x)^2 + 16*
a^2*cos(3*x)^2 + 36*a^2*cos(2*x)^2 + 16*a^2*cos(x)^2 + a^2*sin(4*x)^2 + 16*a^2*sin(3*x)^2 + 36*a^2*sin(2*x)^2
+ 48*a^2*sin(2*x)*sin(x) + 16*a^2*sin(x)^2 + 8*a^2*cos(x) + a^2 + 2*(4*a^2*cos(3*x) + 6*a^2*cos(2*x) + 4*a^2*c
os(x) + a^2)*cos(4*x) + 8*(6*a^2*cos(2*x) + 4*a^2*cos(x) + a^2)*cos(3*x) + 12*(4*a^2*cos(x) + a^2)*cos(2*x) +
4*(2*a^2*sin(3*x) + 3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(4*x) + 16*(3*a^2*sin(2*x) + 2*a^2*sin(x))*sin(3*x)), x)
+ 2*((12*sqrt(2)*x*cos(3/2*x) - (9*sqrt(2)*x^2 - 8*sqrt(2))*sin(3/2*x))*cos(3*x) + 3*(12*sqrt(2)*x*cos(3/2*x)
- (9*sqrt(2)*x^2 - 8*sqrt(2))*sin(3/2*x))*cos(2*x) + 3*(9*sqrt(2)*x^2*sin(x) + 12*sqrt(2)*x*cos(x) + 4*sqrt(2
)*x)*cos(3/2*x) + 9*(3*sqrt(2)*x^2*cos(3/2*x) + 4*sqrt(2)*x*sin(3/2*x))*sin(2*x) - (9*sqrt(2)*x^2 - 36*sqrt(2)
*x*sin(x) + 3*(9*sqrt(2)*x^2 - 8*sqrt(2))*cos(x) - 8*sqrt(2))*sin(3/2*x))*sin(3*x))*sin(2/3*arctan2(sin(3/2*x)
, cos(3/2*x))) - 4*(sqrt(2)*cos(3*x)^3 + 3*(2*sqrt(2)*cos(2*x) + 2*sqrt(2)*cos(x) + sqrt(2))*cos(3*x)^2 + (sqr
t(2)*cos(3*x) + sqrt(2))*sin(3*x)^2 + 9*sqrt(2)*cos(2*x)^2 + 9*sqrt(2)*cos(x)^2 + 9*sqrt(2)*sin(2*x)^2 + 18*sq
rt(2)*sin(2*x)*sin(x) + 9*sqrt(2)*sin(x)^2 + 3*(3*sqrt(2)*cos(2*x)^2 + 3*sqrt(2)*cos(x)^2 + 3*sqrt(2)*sin(2*x)
^2 + 6*sqrt(2)*sin(2*x)*sin(x) + 3*sqrt(2)*sin(x)^2 + 2*(3*sqrt(2)*cos(x) + 2*sqrt(2))*cos(2*x) + 4*sqrt(2)*co
s(x) + sqrt(2))*cos(3*x) + 6*(3*sqrt(2)*cos(x) + sqrt(2))*cos(2*x) + 6*((sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*co
s(3*x) + sqrt(2)*sin(2*x) + sqrt(2)*sin(x))*sin(3*x) + 6*sqrt(2)*cos(x) + sqrt(2))*sin(1/3*arctan2(sin(3/2*x),
cos(3/2*x))))/((a*cos(3*x)^4 + a*sin(3*x)^4 + 2*(3*a*cos(2*x) + 3*a*cos(x) + 2*a)*cos(3*x)^3 + 6*(a*sin(2*x)
+ a*sin(x))*sin(3*x)^3 + 3*(3*a*cos(2*x)^2 + 3*a*cos(x)^2 + 3*a*sin(2*x)^2 + 6*a*sin(2*x)*sin(x) + 3*a*sin(x)^
2 + 6*(a*cos(x) + a)*cos(2*x) + 6*a*cos(x) + 2*a)*cos(3*x)^2 + 9*a*cos(2*x)^2 + 9*a*cos(x)^2 + 9*(a*cos(3*x)^2
+ 9*a*cos(2*x)^2 + 9*a*cos(x)^2 + a*sin(3*x)^2 + 9*a*sin(2*x)^2 + 18*a*sin(2*x)*sin(x) + 9*a*sin(x)^2 + 2*(3*
a*cos(2*x) + 3*a*cos(x) + a)*cos(3*x) + 6*(3*a*cos(x) + a)*cos(2*x) + 6*a*cos(x) + 6*(a*sin(2*x) + a*sin(x))*s
in(3*x) + a)*cos(4/3*arctan2(sin(3/2*x), cos(3/2*x)))^2 + 9*(a*cos(3*x)^2 + 9*a*cos(2*x)^2 + 9*a*cos(x)^2 + a*
sin(3*x)^2 + 9*a*sin(2*x)^2 + 18*a*sin(2*x)*sin(x) + 9*a*sin(x)^2 + 2*(3*a*cos(2*x) + 3*a*cos(x) + a)*cos(3*x)
+ 6*(3*a*cos(x) + a)*cos(2*x) + 6*a*cos(x) + 6*(a*sin(2*x) + a*sin(x))*sin(3*x) + a)*cos(2/3*arctan2(sin(3/2*
x), cos(3/2*x)))^2 + (2*a*cos(3*x)^2 + 9*a*cos(2*x)^2 + 9*a*cos(x)^2 + 9*a*sin(2*x)^2 + 18*a*sin(2*x)*sin(x) +
9*a*sin(x)^2 + 2*(3*a*cos(2*x) + 3*a*cos(x) + 2*a)*cos(3*x) + 6*(3*a*cos(x) + a)*cos(2*x) + 6*a*cos(x) + 2*a)
*sin(3*x)^2 + 9*a*sin(2*x)^2 + 18*a*sin(2*x)*sin(x) + 9*a*sin(x)^2 + 9*(a*cos(3*x)^2 + 9*a*cos(2*x)^2 + 9*a*co
s(x)^2 + a*sin(3*x)^2 + 9*a*sin(2*x)^2 + 18*a*sin(2*x)*sin(x) + 9*a*sin(x)^2 + 2*(3*a*cos(2*x) + 3*a*cos(x) +
a)*cos(3*x) + 6*(3*a*cos(x) + a)*cos(2*x) + 6*a*cos(x) + 6*(a*sin(2*x) + a*sin(x))*sin(3*x) + a)*sin(4/3*arcta
n2(sin(3/2*x), cos(3/2*x)))^2 + 9*(a*cos(3*x)^2 + 9*a*cos(2*x)^2 + 9*a*cos(x)^2 + a*sin(3*x)^2 + 9*a*sin(2*x)^
2 + 18*a*sin(2*x)*sin(x) + 9*a*sin(x)^2 + 2*(3*a*cos(2*x) + 3*a*cos(x) + a)*cos(3*x) + 6*(3*a*cos(x) + a)*cos(
2*x) + 6*a*cos(x) + 6*(a*sin(2*x) + a*sin(x))*sin(3*x) + a)*sin(2/3*arctan2(sin(3/2*x), cos(3/2*x)))^2 + 2*(9*
a*cos(2*x)^2 + 9*a*cos(x)^2 + 9*a*sin(2*x)^2 + 18*a*sin(2*x)*sin(x) + 9*a*sin(x)^2 + 9*(2*a*cos(x) + a)*cos(2*
x) + 9*a*cos(x) + 2*a)*cos(3*x) + 6*(3*a*cos(x) + a)*cos(2*x) + 6*a*cos(x) + 6*(a*cos(3*x)^3 + 3*(2*a*cos(2*x)
+ 2*a*cos(x) + a)*cos(3*x)^2 + 9*a*cos(2*x)^2 + 9*a*cos(x)^2 + (a*cos(3*x) + a)*sin(3*x)^2 + 9*a*sin(2*x)^2 +
18*a*sin(2*x)*sin(x) + 9*a*sin(x)^2 + 3*(3*a*cos(2*x)^2 + 3*a*cos(x)^2 + 3*a*sin(2*x)^2 + 6*a*sin(2*x)*sin(x)
+ 3*a*sin(x)^2 + 2*(3*a*cos(x) + 2*a)*cos(2*x) + 4*a*cos(x) + a)*cos(3*x) + 6*(3*a*cos(x) + a)*cos(2*x) + 6*a
*cos(x) + 3*(a*cos(3*x)^2 + 9*a*cos(2*x)^2 + 9*a*cos(x)^2 + a*sin(3*x)^2 + 9*a*sin(2*x)^2 + 18*a*sin(2*x)*sin(
x) + 9*a*sin(x)^2 + 2*(3*a*cos(2*x) + 3*a*cos(x) + a)*cos(3*x) + 6*(3*a*cos(x) + a)*cos(2*x) + 6*a*cos(x) + 6*
(a*sin(2*x) + a*sin(x))*sin(3*x) + a)*cos(2/3*arctan2(sin(3/2*x), cos(3/2*x))) + 6*((a*sin(2*x) + a*sin(x))*co
s(3*x) + a*sin(2*x) + a*sin(x))*sin(3*x) + a)*cos(4/3*arctan2(sin(3/2*x), cos(3/2*x))) + 6*(a*cos(3*x)^3 + 3*(
2*a*cos(2*x) + 2*a*cos(x) + a)*cos(3*x)^2 + 9*a*cos(2*x)^2 + 9*a*cos(x)^2 + (a*cos(3*x) + a)*sin(3*x)^2 + 9*a*
sin(2*x)^2 + 18*a*sin(2*x)*sin(x) + 9*a*sin(x)^2 + 3*(3*a*cos(2*x)^2 + 3*a*cos(x)^2 + 3*a*sin(2*x)^2 + 6*a*sin
(2*x)*sin(x) + 3*a*sin(x)^2 + 2*(3*a*cos(x) + 2*a)*cos(2*x) + 4*a*cos(x) + a)*cos(3*x) + 6*(3*a*cos(x) + a)*co
s(2*x) + 6*a*cos(x) + 6*((a*sin(2*x) + a*sin(x))*cos(3*x) + a*sin(2*x) + a*sin(x))*sin(3*x) + a)*cos(2/3*arcta
n2(sin(3/2*x), cos(3/2*x))) + 6*((a*sin(2*x) + a*sin(x))*cos(3*x)^2 + 2*(a*sin(2*x) + a*sin(x))*cos(3*x) + a*s
in(2*x) + a*sin(x))*sin(3*x) + 6*(a*sin(3*x)^3 + 6*(a*sin(2*x) + a*sin(x))*sin(3*x)^2 + (a*cos(3*x)^2 + 9*a*co
s(2*x)^2 + 9*a*cos(x)^2 + 9*a*sin(2*x)^2 + 18*a*sin(2*x)*sin(x) + 9*a*sin(x)^2 + 2*(3*a*cos(2*x) + 3*a*cos(x)
+ a)*cos(3*x) + 6*(3*a*cos(x) + a)*cos(2*x) + 6*a*cos(x) + a)*sin(3*x) + 3*(a*cos(3*x)^2 + 9*a*cos(2*x)^2 + 9*
a*cos(x)^2 + a*sin(3*x)^2 + 9*a*sin(2*x)^2 + 18*a*sin(2*x)*sin(x) + 9*a*sin(x)^2 + 2*(3*a*cos(2*x) + 3*a*cos(x
) + a)*cos(3*x) + 6*(3*a*cos(x) + a)*cos(2*x) + 6*a*cos(x) + 6*(a*sin(2*x) + a*sin(x))*sin(3*x) + a)*sin(2/3*a
rctan2(sin(3/2*x), cos(3/2*x))))*sin(4/3*arctan2(sin(3/2*x), cos(3/2*x))) + 6*(a*sin(3*x)^3 + 6*(a*sin(2*x) +
a*sin(x))*sin(3*x)^2 + (a*cos(3*x)^2 + 9*a*cos(2*x)^2 + 9*a*cos(x)^2 + 9*a*sin(2*x)^2 + 18*a*sin(2*x)*sin(x) +
9*a*sin(x)^2 + 2*(3*a*cos(2*x) + 3*a*cos(x) + a)*cos(3*x) + 6*(3*a*cos(x) + a)*cos(2*x) + 6*a*cos(x) + a)*sin
(3*x))*sin(2/3*arctan2(sin(3/2*x), cos(3/2*x))) + a)*sqrt(a))

Giac [F]

\[ \int \frac {x^2}{(a+a \cos (x))^{3/2}} \, dx=\int { \frac {x^{2}}{{\left (a \cos \left (x\right ) + a\right )}^{\frac {3}{2}}} \,d x } \]

[In]

integrate(x^2/(a+a*cos(x))^(3/2),x, algorithm="giac")

[Out]

integrate(x^2/(a*cos(x) + a)^(3/2), x)

Mupad [F(-1)]

Timed out. \[ \int \frac {x^2}{(a+a \cos (x))^{3/2}} \, dx=\int \frac {x^2}{{\left (a+a\,\cos \left (x\right )\right )}^{3/2}} \,d x \]

[In]

int(x^2/(a + a*cos(x))^(3/2),x)

[Out]

int(x^2/(a + a*cos(x))^(3/2), x)